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Abstract. We introduce a formalism for handling general spaces of hierarchical tilings, a 
category that includes substitution tilings, Bratteli-Vershik systems, S-adic transformations, 
and multi-dimensional cut-and-stack transformations. We explore ergodic, spectral and 

£N| topological properties of these spaces. We show that familiar properties of substitution 

tilings carry over under appropriate assumptions, and give counter-examples where these 
assumptions arc not met. For instance, we exhibit a minimal tiling space that is not uniquely 
ergodic, with one ergodic measure having pure point spectrum and another ergodic measure 

Q\ having mixed spectrum. We also exhibit a 2-dimensional tiling space that has pure point 

measure-theoretic spectrum but is topologically weakly mixing. 
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1. Introduction 

Hierarchical structures are ubiquitous in the real world. Typically there are a finite num- 
ber of levels, ranging from the tiny (say, subatomic particles) to the huge (say, clusters of 
galaxies). In many cases the smallest level is so small that it makes sense to extrapolate 
mathematically to infinitely small hierarchical structures - fractals. In this paper we con- 
sider the complementary situation where the smallest scale may not be small, but the largest 
scale is so large that it makes sense to extrapolate to infinite size. 

There is an extensive literature devoted to expanding hierarchies, dating back to the 1800s 
[47], with applications to dynamics dating back to the early 1900s [39]. Most of the aperiodic 
sets of tiles that were discovered over the years, from Berger [11] to Robinson [52] to Penrose 

] to Goodman-Strauss [31] and others, used hierarchy as means of proving aperiodicity. 
Tiles group into clusters that group into larger clusters, etc., so that the resulting patterns 
exhibit structure at arbitrarily large length scales and cannot be periodic. 

In most of the literature, it is assumed that the hierarchies have essentially the same 
structure at each level, so that the system can be described by a single substitution map. 
Indeed, there has been tremendous progress on substitution sequences, substitution subshifts, 
and substitution tilings. However, there is much to be said about hierarchical systems where 
the structure is not necessarily repeated at each level. 
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The idea of studying general hierarchical systems can be seen in the cut-and-stack formal- 
ism of ergodic theory. The first example of Chacon [ ], which exhibited a weakly mixing 
system that was not strongly mixing, is a fusion of the sort discussed in this paper. Over 
the years the technique has been used to construct many interesting examples, and it has 
been shown [5] that all interval exchange transformations, and indeed all aperiodic mea- 
sure preserving transformations, can be obtained by cutting and stacking. Cutting and 
stacking has been generalized to higher dimensions for Z d actions [53, 35], for Reactions 
on rectangular domains [17], and for general locally compact second countable groups [18] 
and amenable groups. Progress has recently been made on nonstationary Bratteli-Vershik 
systems [24, 23, 12], most of which can be viewed as a discrete 1-dimensional version of the 
fusion tilings described in this paper [10]. 

This paper provides a framework for studying the ergodic theory and topology of hier- 
archical tilings. Our formalism encompasses, among other things, substitution tilings and 
substitution subshifts, cut-and-stack transformations, S-adic transformations [ ], and sta- 
tionary and non-stationary Bratteli-Vershik systems [24, 12]. 

Taken to extremes, our formalism can be made too general. Without simplifying as- 
sumptions, essentially any tiling space can be viewed as a fusion, and almost any sort of 
dynamical behavior is possible. For instance, Jewett [34] and Krieger [, ] showed that any 
ergodic measurable automorphism of a non-atomic Lebesgue space system can be realized 
topologically as a uniquely ergodic map on a Cantor set; in most cases these can be viewed 
as subshifts, and hence as fusion tiling spaces. Downarowicz [21] showed that there exist 
Cantor dynamical systems whose invariant measures match an arbitrary Choquet simplex. 

In this paper we identify appropriate hypotheses that preserve the essential properties of 
substitutions while applying to more general systems. Certain properties, like minimality 
or unique ergodicity, hold under very general conditions. Others, like finitely generated 
(rational Cech) cohomology or pure point spectrum or (on the other extreme) topological 
weak mixing, require stronger assumptions. 

In addition, we develop a number of examples that show how these properties can be lost 
when the assumptions are too weak. We hope that these examples will help to classify fusion 
tilings, and to better organize our understanding of tilings in general. 

Some of our proofs are quite simple, yet determining how to apply the techniques of 
substitution systems to fusions is far from trivial. The key tools for studying substitution 
systems are Perron-Frobenius theory and the existence of a self-map that can be iterated 
arbitrarily many times. Neither of these work for general fusions. The new methods devised 
in this paper provide us with more insight into how properties of tiling spaces are related to 
properties of tilings. Some properties of a hierarchical tiling space are directly related to the 
geometry of the individual tiles. Others come from the details of how the tiles are assembled 
into bigger and bigger clusters. Still others can be deduced from coarser numerical data, 
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such as from the matrices that count how many of each kind of tile appear in each kind 
of cluster. Because the hierarchy in fusion rules is less rigid than that of their substitutive 
counterparts, combinatorics, geometry, algebra, and topology can have effects that need to be 
teased apart. The challenge is to understand which properties come from which information, 
and to organize that information effectively. 

Acknowledgments. We thank Mike Boyle, Lewis Bowen, Kariane Calta, Amos Nevo, 
E. Arthur Robinson, Jr. and Boris Solomyak for helpful discussions. The work of L.S. is 
partially supported by NSF grants DMS-0701055 and DMS-1101326. 

2. Definitions 

In this work a prototile is a labelled, closed topological disk in M. d . The label, which can 
be thought of as a color or a marking, is necessary when we wish to distinguish between 
prototiles that are geometrically similar. In general we assume that we have a finite set V of 
prototiles to use as building blocks for our tilings. (This assumption is useful but not entirely 
necessary. In a separate work [: ] we consider tilings built from an infinite but compact set 
P.) We also assume that we have fixed a closed subgroup G of the Euclidean group E(d) 
that contains a full rank lattice of translations; this group G will be used to construct our 
tiles, patches, and tilings and can also serve as the group action of our dynamical system. 
(The two standard translation subgroups that appear in tiling theory are Z d and R d .) It is 
possible to act on a prototile by an isometry in G by applying the isometry to the closed set 
defining the prototile and carrying the labelling information along unchanged. A prototile 
which has been so moved is called a tile. We will abuse notation by denoting the application 
of an isometry g E G to a prototile p as g(p); when the isometry is translation by t? G M d we 
denote the translated tile by p + v. A V -patch (or patch, for short) of tiles is a connected, 
finite union of tiles that only overlap on their boundaries; the support of the patch is the 
closed set in W 1 that it covers. Two tiles or patches are considered equivalent or copies of one 
another if there is an element of G taking one to the other. A tiling T of M. d is a collection 
of tiles that completely cover R d and overlap only on their boundaries. 

A tiling is said to have finite local complexity (FLC) with respect to the group G if it 
contains only finitely many connected two-tile patches up to motions from G. Most of 
the literature on tiling dynamical systems uses finite local complexity as a key assumption. 
This work in this paper is limited to FLC fusion tilings. Fusion tilings with infinite local 
complexity (LLC) will be considered in [29]. 

2.1. Fusion tilings. Given two P-patches Pi and P2 and two isometries g\ and g<i in G, if the 
patches g\{P\) and #2(^2) overlap only on their boundaries, and if the union #i(Pi) U #2(^2) 
forms a P-patch, we call that union the fusion of Pi to P2 via g\ and #2- When we do not 
wish to specify the isometries we may call it a fusion of Pi to P2. Notice that there will be 
many ways to fuse two patches together and that we may attempt to fuse any finite number 
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of patches together. We may even fuse a patch to copies of itself. Patch fusion is simply a 
version of concatenation for geometric objects. 

The idea behind a "fusion rule" is an analogy to an atomic model: we have atoms, and 
those atoms group themselves into molecules, which group together into larger and larger 
structures. In this analogy we think of prototiles as atoms and patches as molecules. Let 
Po = P be our prototile set, our "atoms". The first set of "molecules" they form will be 
defined as a set of finite P-patches Pi, with notation V\ = {P\(l), -Pi (2), P\{jx)}. Next 
we construct the structures made by these "molecules": the set V2 will be a set of finite 
patches that are fusions of the patches in V\. That is P 2 = {-fMl), -^2(2), P2O2)} is a 
set of patches, each of which is a fusion of patches from V\. While the elements of V2 are 
technically P-patches, we can also think of them as Pi-patches by considering the elements 
of Pi as prototiles. We continue in this fashion, constructing P 3 as a set of patches that 
are fusions of patches from V2 and in general constructing P n as a set of patches which are 
fusions of elements of P n -i- The elements of P n are called n-fusion supertiles or n-supertiles, 
for short. 1 We collect them together into an atlas of patches we call our fusion rule: 

TZ = {V n ,n GN} = {P n {j) I n e N and 1 < j < j n } . 

A patch is admitted by TZ if a copy of it can be found inside some supertile P n (j) for some 
n and j. A tiling T of lR d is said to be a fusion tiling with fusion rule TZ if every patch of tiles 
contained in T is admitted by TZ. We denote by X-r the set of all P.-fusion tilings. Given 
a fusion rule, we can obtain another fusion rule TZ' with j' n = j n+1 and P' n {j) = P n+ i(j). 
We simply ignore the lowest level and treat the 1-fusion supertiles as our basic tiles. The 
resulting tiling space is denoted X^. Likewise, X\ is the space of tilings obtained from TZ 
in which the /c-fusion supertiles are considered the smallest building blocks. 

Standing assumption (for this entire paper): If none of the supertiles in TZ have 
inner radii approaching infinity then X-ji will be empty, so for that reason we restrict our 
attention to fusion rules that have nontrivial tiling spaces. 

When d — 1 and G = Z, with all tiles having unit length, fusion tilings correspond to 
Bratteli-Vershik systems, modulo complications having to do with edge sequences that have 
no predecessors or no successors. See [10] for more about the correspondence. (In addition 
to subshifts, Bratteli-Vershik systems can model non-expansive maps on Cantor sets; these 
can also be viewed as 1-dimensional fusion tilings, albeit with infinitely many tile types [29].) 

Example 2.1. The Chacon transformation. In [14] there is an early example of a transfor- 
mation that is weakly mixing but not strongly mixing. The original cutting-and-stacking 
construction is a self-map on an interval; the stacking portion can be seen as a sort of fusion. 

^If we wish, we can also add labels to the supertiles, so that the information carried in an n-supertile is 
more than just its composition as a patch in a tiling. This generalization is useful for collaring constructions, 
as in Section 5. 
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However for the purposes of an immediate example we use the fact that the Chacon space 
can be viewed symbolically using the substitution rule 

a — > aaba b — >■ b, 

which can be iterated by substituting each letter and concatenating the blocks. If we begin 
with an a we have: 

a — y aaba — > aaba aaba b aaba —>■... 

In order to make a Chacon tiling of K. we only need to assign closed intervals to the symbols 
a and b and place them on the line according to the symbols in a Chacon sequence. 

We can view a Chacon tiling of R as a fusion tiling as follows. Consider l a and Ij, to be two 
positive numbers and let a denote a prototile with support [0, l a ] and b denote a prototile 
with support [0,Z&]. (If l a = h then we use the symbols a and b as labels to tell the tiles 
apart) . We define Pi (a) = a U (a + l a ) U {b + 2l a ) U (a + 2l a + l b ) and Pi (b) = b. The length of 
Pi (a) is 3l a + lb- To make P2(a) we simply fuse three copies of Pi (a) and one copy of Pi(b) 
together in the correct order, and of course P2Q)) = b still. The length of the new a supertile 
is three times that of the previous a supertile plus the length of b. We continue recursively 
to construct all of the n-fusion supertiles. 

2.1.1. Transition matrices and the subdivision map. Given a fusion rule 1Z there is a family 
of transition matrices that keep track of the number and type of (n — l)-fusion supertiles 
that combine to make the n-supertiles. The transition matrix for level n, denoted M n _i >n , 
has entries M n _i jn (fc, /) = the number of (n — l)-supertiles of type k, that is, equivalent to 
P n -i(k), in the n-supertile of type /, P n (i). If there is more than one fusion of P ra _i-supertiles 
that can make P n (l), we fix a preferred one to be used in this and all other computations. 
For levels n < N e N, we likewise define the transition matrix from n- to iV-supertiles as 
M^jsr = M n ^ n+ iM n+ i t n + 2 • • • Mjv-i,iv- The entry of M n ^ is the number of n-supertiles 
of type % in the iV-supertile of type j. Another way to think about this is to imagine a 
"population vector" v G 2P N of a patch of iV-supertiles: the entries represent the number of 
iV-fusion supertiles of each type appearing in the patch. Then M N _ ±>N v gives the population 
of this patch in terms of (iV — l)-supertiles, Mn^2,n-iM^_i^v gives the population in terms 
of (N — 2)-supertiles, and M n ^v gives the population of this patch in terms of n-supertiles. 

Any self-afnne substitution tiling, in any dimension, can be viewed as a fusion tiling. An 
n-supertile is what we get by applying the substitution n times to an ordinary tile, and 
can be decomposed into (n — l)-supertiles according to the pattern of the substitution. For 
such tilings, the matrix M n ^ is just the (N — n)th power of the usual substitution matrix. 
However, there is an important difference in perspective between substitutions and fusions. 

A substitution can be viewed as a map from a tiling space to itself, in which all tiles are 
enlarged and then broken into smaller pieces. This map can be repeated indefinitely. In 
a fusion tiling, we can likewise break each n-fusion tile into level (n — l)-supertiles using 
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the subdivision map a n , which is a map from X'^ to X^ . Unlike the substitution map 
for self-affine tilings, it cannot go from X-ji to itself, and this map cannot be repeated more 
than n times. Once you are down to the atomic level (i.e., ordinary tiles), you cannot 
subdivide further! The proofs of theorems about substitution tilings often involve taking 
an arbitrary tiling and applying a substitution, or sometimes its inverse, enough times to 
achieve a desirable result. For general fusion tilings, this line of reasoning usually does not 
work. 

2.1.2. Induced fusions. Let {N(n)}^ =l be an increasing sequence of positive integers. The 
induced fusion on N(n) levels, lZ md , is obtained from a given fusion 1Z by composing the 
fusions for levels N(n) + 1, N(n + 1) into one step. In this case the supertiles of TZ md are 
given by V™ d = Pjv(n), where the iV(n)-supertiles are seen as fusions of N(n — l)-supertiles. 
The transition matrices for 'R, ind are given by M™ n+1 = M N ( n ^ N ( n+ iy 

2.1.3. All FLC tilings are fusion tilings. It is possible to view any tiling T of M. d from a given 
prototile set Vq as a fusion tiling, as long as it has finite local complexity. Let the set V n 
consist of all connected patches containing n tiles or less. By finite local complexity this is 
a finite set. Each element of V n is either an element of V n -\ or is the fusion of two elements 
of V n -\. (The fact that these fusions typically are not unique does not matter). 

2.2. Common assumptions. The previous section shows that the category of fusion tilings 
is extremely general. To prove meaningful results, we have to impose additional conditions 
on our fusion rules. We collect several of them into this section. 

2.2.1. Prototile- and transition-regularity. These are the cases that are most similar to the 
usual definitions of symbolic and tiling substitution. When the number of supertiles at 
each level is constant, we can associate each n-supertile to a specific prototile, regardless 
of whether there is a geometric connection between the two. When we do this we call the 
fusion rule prototile-regular and rewrite it as: 

Tl = {P n (p) | n E N and p E V } . 

If the number j n of supertiles at the nth level of a fusion rule 1Z has J = liminf j n for 
some finite J, then the fusion rule is equivalent to a prototile-regular fusion rule by inducing 
on the levels for which j n = J. The price we pay for taking such an induced fusion is that 
the transition matrices can become wildly unbounded. 

In the special case where the number of supertiles at each level is a fixed constant J, if the 
transition matrices are all equal to a single matrix we call the fusion rule a transition-regular 
fusion rule. Being transition-regular is considerably stronger that being prototile-regular. 
All substitution sequences and self-affine tilings as defined in, for instance, [50, 57] are 
transition-regular, but not every transition-regular fusion tiling comes from a substitution. 
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The combinatorics and geometry of how the (n — l)-supertiles join to form n-supertiles can 
change from level to level. 

Example 2.2. A fusion that is transition-regular but not a substitution. Consider a 1- 
dimensional fusion rule with transition matrix (\ \) in which P n (a) is always given by the 
word P n _i(a)P n _i(a)P n _i(6), and in which P n (b) is given by P n _i(&)P n _i(6)P n _i(a) if n is 
prime, and is given by P n _i(a)P rt _i(6)P„_i(6) if n is composite. 

Remarks. 

(1) Pseudo-self-similar (or self-amne) tilings, such as the Penrose tiling with kites and 
darts, are also transition-regular fusion tilings. In many cases these are asymptoti- 
cally self-similar, and this asymptotic structure was used [27, 49] to show that such 
tilings were topologically equivalent to self-similar tilings with fractal boundaries. 

(2) In the correspondence between one-dimensional fusion tilings with G = Z and Bratteli- 
Vershik systems, prototile-regular tilings correspond to finite Bratteli diagrams. The 
finite list of vertices on the nth level of the Bratteli diagram represents the finite set 
of n-supertiles. 

(3) The one- dimensional S-adic substitution sequences of Durand [ !] can be recast as 
fusion tilings, as can the linearly recurrent Delone sets and tower systems in [9, 3]. 
Example 2.2 is 5-adic. 

(4) The "non-constructive" combinatorial substitutions in [25] are exactly the class of 
prototile-regular fusion tilings. 

2.2.2. Primitivity. A fusion rule is said to be primitive if, for each non-negative integer n, 
there exists another integer iV such that every n-fusion supertile is contained in every N- 
supertile. When the fusion rule is transition-regular this is equivalent to some power of the 
transition matrix having strictly positive entries. In general it is equivalent to there existing 
an N for each n such that M n ^ has all positive entries. A fusion rule is called strongly 
primitive if for every n > 1, each (n + l)-supertile contains at least one copy of every n- 
supertile. That is, all of the transition matrices M„ in+ i have strictly positive entries. Any 
primitive fusion rule is equivalent to a strongly primitive one by inducing on enough levels. 

Primitivity is one of the most common assumptions used in the literature on substitution 
sequences and tilings. It allows for Perron-Frobenius theory to be applied to the systems to 
determine natural frequencies, volumes, and expansion rates. We will adapt this analysis to 
the fusion situation in Section 3. 

2.2.3. Recognizability. A fusion rule 71 is said to be recognizable if, for each n, the subdivision 
map a n from A^ to A^T 1 is a homeomorphism. If so, then every tiling in A^ can be 
unambiguously expressed as a tiling with n-supertiles for every n. The uniform continuity of 
the inverse subdivision maps then implies that there exists a family of recognizability radii 
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r n (n = 1, 2, . . .), such that, whenever two tilings in X-ji have the same patch of radius r n 
around a point v G M d , then the n-supertiles intersecting v in those two tilings are identical. 

For substitution sequences and tilings, recognizability is closely related to non-periodicity 
[40, 58]. Recognizability implies that none of the tilings are periodic. Conversely, if G 
consists only of translations [58], or if G contains a set of rotation about the origin with no 
invariant subspaces then the absence of periodic tilings in X-ji implies recognizability [. ]. 
However, it is easy to construct fusion rules that are nonperiodic but not recognizable. For 
instance, the Fibonacci tiling can be generated either from the substitution a — > ab, b — > a 
or from the substitution a — > ba, b — > a. By including both sets of supertiles in our fusion 
rule, we obtain a description of the non-periodic Fibonacci tiling space in which each tiling 
has at least two (actually more) decompositions into n-supertiles for n > 0. 

We now show that fusion tiling spaces are topological factors of recognizable fusion tiling 
spaces using a construction inspired by the work of Robinson [52] and Mozes [41]. 

Example 2.3. Constructing a recognizable extension. Let TZo be a 1-dimensional fusion rule 
on the letters a and b, each of which is viewed as a tile of length 1. If we let the n-supertiles 
be all possible sequences of as and 6s of length 5 n , then the space Xn is just the space of 
all bi-infinite tilings by a's and 6's and TZq is clearly not recognizable. 

Now let 1Z be a 1-dimensional fusion with four letters, a\, a 2 , b\ and b 2 - We call a\ and 
b\ "type 1", and write x\ to mean either a\ or b\. Likewise X2 means either 02 or 62- The 
1-supertiles are all 5-letter words of the general form x 2 XiX±XiXi (where each Xi denotes a 
separate choice of a, or bi) or x 2 XiX 2 x 2 Xi. We will use s{ are shorthand for supertiles of the 
first type and si for the second. Note that each supertile begins with an isolated x 2 , and 
that isolated x 2 s appear only at the beginning of supertiles. This makes the map from X\ 
to X-ji invertible. 

We repeat the coding at higher levels. Second-order supertiles can either take the form 
s\ = sls{s{s{s{ or si = sls\slsls\, and generally (n + l)-supertiles can either take the form 
s'1 +1 = s^s^s^SiS^ or S2 +1 = s^s^s^s^s™. By the same reasoning, all decomposition maps 
are invertible, and 71 is recognizable. Finally, the factor map Xn — > Xn just erases the 
subscripts on all of the letters. 

The details of the construction will be different for different examples and can get compli- 
cated if the supertiles have wild shapes or combinatorics, but the basic idea is universal. Pick 
sufficiently many copies of your original tile set. Use some of the labels within a first-order 
supertile to indicate which tiles are in the supertile, and the rest to give the first-order super- 
tiles labels. Use some of those first-order labels to define the boundaries of the second-order 
supertiles, and the rest to label the second-order supertiles. By continuing the process ad 
infinitum we obtain a recognizable fusion tiling space that factors onto the original. How 
close to one-to-one this factor map becomes an important question. 
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2.2.4. Van Hove sequences and fusion rules. A van Hove sequence {An} of subsets of M. d 
consists of sets whose boundaries are increasingly trivial relative to their interiors in a precise 
sense. In many cases it will be convenient to consider only fusion rules where the supertiles 
share this property. The use of van Hove sequences, which for M. d is equivalent to F0lner 
sequences, is adopted from statistical mechanics. We follow the notation of [57] here and 
define, for any set A G R and r > 0: 

A +r = {xeM d : dist(a;, A) < r}, 

where "dist" denotes Euclidean distance. A sequence of sets {A n } of sets in R d is called a 
van Hove sequence if for any r > 

ii m vd 'y ) =°. 

ra->oo \/0l{A n ) 

where dA is the boundary of A and Vol is Euclidean volume. 

Given a fusion rule 1Z, we may make a sequence of sets in M. d by taking one n-supertile 
for each n and calling its support A n . We say 1Z is a van Hove fusion rule if every such 
sequence is a van Hove sequence. Equivalently, a fusion rule is van Hove if for each e > 
and each r > there exists an integer n such that each n-supertile A, with n > n , has 
Vol{dA) +r < eVol(A). 



2.3. Notational conventions. Entries of vectors and matrices are indicated as arguments, 
while subscripts are used to distinguish between different vectors and matrices. Thus, 
M 12 (3,4) is the (3,4) entry of the matrix M 12 and 1*5(2) is the second entry of t>5. Vec- 
tors are viewed as columns, so that the product Mv of a matrix and a vector is well-defined. 
Groups are denoted by capital letters, as are subsets of groups, while elements of groups are 
lower case. Collections of patches of tilings are given by calligraphic letters V, TZ, etc, and in 
particular our fusion rules are so denoted. Tilings are bold face. Elements of physical space 
M. d are marked with arrows, and the dot product is reserved for this setting. 



3. Dynamics of fusion tilings 

Let G = Gt x G r , where G t is the translation subgroup and G r is the point group G/G t - 
By assumption G contains a full rank lattice of translations, and G r is a closed subgroup of 
0(n). Let Vol be Haar measure 011 G fl a product of Lebesgue measure in the continuous 
directions of G t and counting measure in the discrete directions, and let Ao be normalized 
Haar measure on G r . Let A be a measure on G with X(U t x U r ) = Vol(U t )\o(U r ) for every 
pair of measurable sets U t G G t , U r G G r . We assume that we have a metric on G whose 
restriction to Gt is Euclidean distance and whose restriction to G r is bounded by 1. 
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3.1. Tiling metric topology and dynamical system. Let X v be the set of all possible 
tilings using some fixed prototile set V and some group G of isometries. (That is, a point 
in this space is an entire tiling of W 1 .) We turn X v into a metric space using the so-called 
"big ball" metric using the metric on G as follows. Two tilings Ti and T 2 are e-close if there 
exist group elements g\ and #2, each of size less than or equal to e, such that <?i(Ti) and 
92(^2) exactly agree on the ball of radius 1/e around the origin. 

This metric is not G-invariant, as it gives greatest weight to points close to the origin, 
but the resulting topology is G- invariant. A sequence of tilings converges to a tiling T if 
there exists a sequence of group elements converging to the identity, such that for every 
compact subset K of Euclidean space, the tilings <7i(Tj) eventually agree with T on K. 

Definition 3.1. Let G' C G contain G t and let X be a closed, G' -invariant subset of X-p. 
A tiling dynamical system (X,G') is the set X together with the action of G' on X. 

It is usually assumed in the tiling literature that G' — Gt- This can be assumed without 
loss of generality when G r is a finite group simply by making extra copies of each prototile, 
one for each element of G r . The situation is more complicated in cases such as the pinwheel 
tiling, where G r is infinite. 

3.2. Minimality. A topological dynamical system (X, G') is minimal if X is the orbit clo- 
sure of any of its elements. 

Proposition 3.2. If the fusion rule 1Z is primitive, then the dynamical system (X n ,G) is 
minimal. 

Proof. Let T in X-ji be any fixed tiling. We will show that given T' G X-ji and e > there 
is a group element g such that which d(g(T),T') < e. Denote by [T'] r the patch of tiles in 
T' that intersect the ball of radius r centered at the origin. By definition we know that any 
such patch is admissible by 1Z and so there is an n G N and a i G {1, j n } for which [T']i/ e 
is a subpatch of P n (i). 

On the other hand, primitivity means that there is an N such that every A^-supertile 
contains a copy of P n (i). Since T is a union of A^-supertiles, it contains many copies of 
P n (i). Pick g to bring some particular copy of P n (i) to the origin in agreement with [T']i/ e . 
Since T' and g(T) are identical on the ball of radius 1/e about the origin, d(g(T), T') < e. □ 

Remarks. 

(1) It is not necessarily true that (X-n,G t ), i.e. the dynamical system with only transla- 
tions acting, is minimal. Consider any fusion rule having only finitely many relative 
orientations of the prototiles, but which for some reason we took G to be the full 
Euclidean group. In this case (X, G) would be minimal but (Xn, Gt) would not. No 
tiling could approximate an irrational rotation of itself. 
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(2) On the other hand, the pinwheel tilings [48] provide an example where G is the full 
Euclidean group but (X, G t ) is minimal. 

(3) Primitivity is sufficient but not necessary for minimality. In particular, the Chacon 
transformation is not primitive, but is minimal. For each n there does not exist an 
N for which P/v(6) contains P n {o). However, there does exists a radius r n such that 
every ball of radius n contains at least one P n (a) and at least one P n (b), so the patch 
[T']i/ e can be found in every T. 

3.3. Invariant measures in general tiling dynamical systems. We begin our treatment 
of the ergodic theory of fusion tilings with a discussion of how invariant measures work for a 
general FLC tiling dynamical systems X, with a focus on patch frequency. For convenience, 
we assume that our action is by G = M. d only. 

Let P be any patch of tiles containing the origin. Let U be a measurable subset of IR d , 
let Xpjj be the cylinder set of all tilings in X that contain P — v for some v e U, and let 
Xp,u to be the indicator function of this cylinder set. The sets Xp t u, plus translates of these 
sets, generate our a-algebra of measurable sets in X. Let \x be an invariant measure on X. 
If U is sufficiently small that, for every tiling Tel, there is at most one v £ U for which 
P — v C T, then ii(Xp t u) must be proportional to the volume of U and we define 

(3-1) freq,(P) = ^l—^Xpfj), 

a quantity that is independent of U. 

For any A C M. d we denote the number of times an equivalent copy of P appears in T, 
completely contained in the set A, as #(P in A R T). As a special case, if P' is another 
patch (usually some supertile), we denote by #(P in P') the number of equivalent copies of 
P completely contained in P' . Next we pick a specific Uq that is a small ball centered at the 
origin and define the function 

(3.2) /p(T) = y^X^m 

This is a smeared 5-function that counts the appearances of P. / /p(T — v)dv is essentially 

J A 

#(P inAnT), with small corrections for patches that come within the diameter of Uq of 

r i 

the boundary of A. Note that / fp(T)dfi = . n(X PUo ) = freq^P). 

Jx Vol[U ) 

We use the following version of the pointwise ergodic theorem: 

Theorem 3.3. Let (X, IR d ) be a tiling dynamical system with invariant Borel probability 
measure \i. Let {A m } be a sequence of balls centered at the origin, with radius going to 
infinity, and let P be any finite patch. Then for ^-almost every tiling T the limit 

(3-3) lim — / f P (T - v)d\(v) = f P (T) 

m->oo V0l{A m ) J Am 
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exists. Furthermore, / /p(T)d/x = / fp( r £)dfi = freq^P). If [i is ergodic, then for 

J x J x 

H-almost every T, /p(T) = freq^(P). 

The quantity /p(T) corresponds to the usual notion of frequency as "number of occur- 
rences per unit computed with an expanding sequence of balls around the origin. 
We will use the term "frequency" for this quantity, and will call freq^(P) an "abstract fre- 
quency". The ergodic theorem says that almost all tilings have well-defined frequencies, and 
that the abstract frequency freq^P), while not necessarily the frequency of P in any specific 
tiling, is the average over all tilings of the frequency of P. Thus, any upper or lower bounds 
on the frequency of P that apply to //-almost every T result in upper or lower bounds on 
the abstract frequency. If every tiling has the same frequency of P, then there is only one 
possible value for the abstract frequency of P, and thus for the measure of any X P)U . Tiling 
spaces where all tilings have the same set of frequencies are uniquely ergodic. 

For an FLC tiling, the set of all patches (up to translation) is countable, and the intersec- 
tion of a countable number of sets of full measure has full measure. As a result, /i-almost 
every tiling T has well-defined frequencies for every patch P. 

Conversely, if a tiling T has well-defined patch frequencies, then we can construct a prob- 
ability measure on X by taking fj,(Xp t u) = f P (T)Vol(U) for small U and extending by 
additivity to larger {7's. The additivity properties of measures follow from the addititiv- 
ity of frequencies. For instance, if a patch P can be extended in two ways, to P' or P", 
then Xp u = X P i tU J J Xpn y. The identity fi(X PiU ) = n(X P i tU ) + fi(X P n jj) follows from 
/p(T) = /p'(T) + fp"(T). Countable additivity is not an issue, since the tiling space is lo- 
cally the product of Euclidean space (where Lebesgue measure has all the desired properties) 
and a Cantor set (where the cr-algebra is based on finite partitions into clopen sets). 

A measure defined in this way may or may not be ergodic. For instance, if T is a one- 
dimensional tiling with the pattern a 00 ^ 00 , with a tiles to the left of the origin and b tiles 
to the right, then the resulting measure on the orbit closure of T is the average of the two 
ergodic measures. 

3.4. Invariant measures and fusion tilings. The possibilities for invariant measures of 
fusion tilings are intimately connected to the asymptotic behavior of the transition matrices 
M n<N as iV — > oo. Our analysis of these matrices takes the place of the standard Perron- 
Frobenius theory used so fruitfully for substitution systems. The results of this subsection 
and the next closely parallel those of [24, 12], the difference being that those papers consider 
discrete systems in 1 dimension, while we consider continuous systems in an arbitrary number 

2 Ergodic theorems are often stated not with balls, but in terms of F0lner or van Hove sequences that 
have special properties, such as being "regular" or "tempered". That generality is useful for computing 
frequencies using different sampling regions, or when considering more complicated groups than K d . For our 
purposes, however, balls are sufficient. 
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of dimensions. For convenience we continue to assume that G = l$L d , and we also assume 
that our fusion rule is van Hove, recognizable, and primitive. 

We define the frequency fp n fj) of a supertile P n (j) in a tiling T to be its frequency as a 
supertile, not as a patch. In other words, /p n m(T) is obtained by viewing T as an element 
of X^, thereby excluding patches that have the same composition as P n {j), but are actually 
proper subsets of another n-supertile or straddle two or more n-supertiles. The abstract 
supertile frequency of P n (j) is obtained by averaging fp n (j) over all tilings. By recognizability, 
each occurrence of a supertile P n {j) is marked by an element of a set of larger patches 
Si. We then have /p„y)(T) = and the abstract supertile frequency of P n (j) is 

i 

i 

Consider a sequence p = {p n } where each p n e M Jn has all nonnegative entries. We say 

jn 

that p is volume-normalized if for all n we have p n {i)V ol{P n {i)) = 1. We say that it has 

i=l 

transition consistency if p n = M u ^Pn whenever n < N. A transition-consistent sequence p 
that is normalized by volume is called a sequence of well-defined supertile frequencies. This 
terminology will be justified by the proof of Theorem 3.4. 

Theorem 3.4. Let 1Z be a recognizable, primitive, van Hove fusion rule. There is a one-to- 
one correspondence between the set of all invariant Borel probability measures on (Xn, M. d ) 
and the set of all sequences of well-defined supertile frequencies with the correspondence that, 
for all patches P, 

jn 

(3.4) freq,(P) = lim V # (P in P n (i)) p n {i) 

i=l 

Proof. Suppose p is an invariant measure. For each n6N and each % e {1,2, j n }, define 
p n (i) to be the abstract supertile frequency of P n (i). For a fixed n, Xn is the union of cylinder 
sets given by which n-supertile is at the origin. Since the measure of X-ji is 1 and the measure 
of each of these cylinder sets is p n (i)Vol(P n (i)), the sequence p is volume-normalized. 

The set of tilings where the origin lies in an n-supertile of type i is the union of disjoint 
sets where the origin lies in a supertile of type i, which in turn sits in an particular way 
in an iV-supertile. There are M n ^(i,j) ways for P n (i) to sit in P^(j). The additivity of 

Jn 

the measure implies that p n {i) = M n ^(i, j)pv(j)- Hence p, gives rise to a sequence of 

j'=i 

well-defined supertile frequencies. 

To see that equation (3.4) applies, let P be any patch of diameter L P . Since the fusion 
rule is van Hove, we can pick an n large enough that the fraction of each n-supertile within 
Lp of the boundary is so small that P patterns appearing in this region can only contribute 
e or less to the frequency of P's in a union of n-supertiles. 



FUSION: A GENERAL FRAMEWORK FOR HIERARCHICAL TILINGS OF R d 



15 



To count the number of P's in a large ball A m around the origin, we must count the number 
of P's in each n-supertile contained in that ball, plus the number of P's that straddle two 
of more n-supertiles, plus the P's in an n-supertile that is only partially in the ball. As 
a fraction of the whole, the third set goes to zero as m — > oo and the second goes to 

jn 

zero as n — > oo. Thus #(P in A m fl T) = ^#(P in P n (i))#(P n (i) in A m )+ boundary 

i=i 

occurrences. Dividing by Vol(A m ) and taking limits, first as m — > oo and then asn-} oo, 
gives the identity 

/ P (T) = lim # ( P in p n(^)) /p n(i )(T) 

i=l 

for all tilings T with well-defined patch frequencies. Integrating this identity over all tilings 
then gives equation (3.4). 

Now suppose that {p n } is a sequence of well-defined supertile frequencies. To establish 
the existence of an invariant measure \i for which {p n } represents the abstract supertile 
frequencies, we simply define freq^(P), and hence the measure of each cylinder set X P>U , by 
equation (3.4). 

To see convergence of the limit on the right hand side, note that, if n < N, the number 
of P in Pjv(j) is the sum of the number of P in each n-supertile in P/v(j), plus a small 
contribution from P's that straddle two or more supertiles. That is, #(P in P/v(j)) w 
^#(P in P n (i))M n , N (i,j), so 

i 

J2*(P ^ Pn(j))pnU) « ]T#(P in P n (i))M n , N (t,j)p N (j) = ^#(P in P n (i))p n (i). 

j i,j i 

As n — > oo the contribution of P's that straddle two n-supertiles goes to zero, so the right 
hand side of (3.4) is a Cauchy sequence. 

The non-negativity of the measure follows from the non- negativity of each p n . The identity 
p(Xn) = 1 follows from volume normalization. Finite additivity follows from the observation 
that, if a patch P can be extended to P' or P" , then #(P in P n {i)) = #(P' in Pn{i)) + 
#(P" in P n (i)), plus a small correction for the situations where P is completely contained 
in P n (i) but P' or P" is not, a correction that does not affect the limit as n — > oo. As 
noted earlier, countable additivity is not an issue for tiling spaces. Thus \i is a well-defined 
measure. 

□ 

3.4.1. Parameterization of invariant measures. In Theorem 3.4 we showed how measures re- 
late to well-defined sequences of supertile frequencies. We now give an explicit parametriza- 
tion of the invariant measures in terms of the transition matrices M n ^, a parametrization 
that we will then use to investigate unique ergodicity. 
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The direction of each column of M n ^ is defined to be the volume-normalized vector 
collinear with it, and we define the direction matrix D n ^ to be the matrix whose columns 
are the directions of the columns of M n ^. That is, 

n ( u\ M n>N (*,k) 

DnN{*,k) = ; ■ . 

ElllMn,N(l,k)V0l(P n (l)) 

Let A n jy be the convex hull of the columns of D n ^, sitting within the hyperplane of 
volume- normalized vectors in W n . Note that the extreme points of A n ^ are columns of 
D n N, but not every column need be an extreme point. Since each column of M n ^ +1 is a 
sum of columns of M n>Nl each column of D n N+ i is a weighted average of the columns of 
D n ,N, so A n> jv+i C A n> jv- Let A n = n^ =n+1 A ni Ar. 

The matrix M n ^ defines an affine map sending An to A n , since if pjy is volume-normalized 
in W N , then so is M n ^PN G W n . We define Aoo to be the inverse limit of the polytopes A n 
under these maps. 

Corollary 3.5. Let (X^,]R d ) be the dynamical system of a recognizable, primitive, van Hove 
fusion rule. The set of all invariant Borel probability measures is parameterized by A^. 

Proof. By Theorem 3.4, we need only show that each element of A M gives rise to a sequence 
{p n } of well-defined supertile frequencies and vice versa. 

By construction, each point in A^ is a sequence of well-defined supertile frequencies, 
since each point in A n is volume-normalized and non-negative, and since the sequence has 
transition consistency. For the converse, suppose that {p n } is a sequence of well-defined 
supertile frequencies. We must show that p n G A n . Since p n = M^jvPjv, Pn is a non-negative 
linear combination of the columns of M n ^, and so is a weighted average of the columns of 
D n N . Thus p n G A n A r. Since this is true for every N, p n G A n . □ 

3.4.2. Measures arising from supertile sequences. In this section we provide a concrete way 
of visualizing certain invariant measures, in particular the ergodic ones. The way to do it is 
by looking at frequencies of patches as they occur in specific sequences of nested supertiles. 

Definition 3.6. Let k = {k n } be a sequence of supertile labels, with k n G {1,2, .. . ,j n }. For 
each n < N , we consider the frequency of each n-supertile P n {i) within P/v(^jv)-' 

Pn,N{i) = M nyN (i,k N )/Vol(P N (k N )). 

We say that k has well-defined supertile frequencies if p n (i) = lim Pn^ii) exists for every 

N— >oo ' 

n and every i G {1, . . . , j n }. 

Note that the vectors p n (i), if they exist, do indeed form a sequence of well-defined supertile 
frequencies. For n < n < N, p n ^ = M n ^ n ip n r^. Taking a limit as N — )■ oo gives p n = 
M n>n ip n /, so the sequence has transition consistency. Likewise, it is easy to check volume 
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normalization. We can therefore associate an invariant measure to every sequence k that 
has well-defined supertile frequencies. 

The purpose of using a sequence k is to visualize a measure. Given such a sequence, 
one can imagine a tiling T where the origin sits inside a k\ 1-supertile, which sits inside 
a k 2 2-supertile, etc. Under mild assumptions, the supertile frequencies /p n m(T) will then 
equal p n (i), and for any patch P, /p(T) will equal freq^P), where \i is constructed from 
the sequence {p n }- The concept of using sequences k to obtain measures applies even to 
non-primitive fusions, as long as the supertile frequencies are well-defined. 

Example 3.7. A minimal fusion rule with two ergodic measures. This is a variation on 
an example found in [23] and illustrates the results of [ ]. Consider a prototile-regular 
1-dimensional fusion rule with two unit length tiles a and b and let G — R. Let P n (a) = 
(-Rn-i(fl)) 10 Pn-i{b) and P n (b) = (P n _i(o)) 10 P n _i(a), so that P\{a) = aaaaaaaaaab and 
Pi{b) = bbbbbbbbbba, etc. M n _ 1)Tl = ( 10 ™ l ^ n ) which has eigenvalues 10 n - 1 and 10™ + 
1. Elementary linear algebra allows us to compute the frequencies as follows. Let a n = 

J^J — ^ , which approaches a limit of just over 0.8 asn-> oo. The fraction of a's in P n (a) 

k=l ~^ 

is (1 + a n )/2 0.9, while the fraction of a's in P n (b) is (1 — a n )/2 « 0.1. 

There are exactly two ergodic measures on this system. A n is an interval for every value 
of n, with endpoints defined by the limits of the first and second columns of D n ^. Likewise, 
Aqo is an interval, whose endpoints fj, a and /i& can be obtained from the supertile sequences 
K = (a, a, a, a, . . .) and k = (6, 6, 6, 6, . . .). The first ergodic measure, fi a , sees the frequencies 
of a's and 6's as measured in the type-a supertiles and thus is a-heavy; the second, /if,, 
reverses the roles of a and b and is b- heavy. The measure w = (fi a + //&)/ '2 is invariant but 
not ergodic; this measure corresponds to Lebesgue measure when the system is seen as a 
cut-and-stack transformation. 

In general, a prototile-regular fusion with j species of tiles can have at most j ergodic 
measures. Of course, there can be fewer, if one or more columns of D n ^ are in the convex 
hull of the others for large N. The following example shows how a sequence k may lead to 
a measure that is not ergodic. 

Example 3.8. A non-ergodic measure coming from a sequence k. Consider the following 
variant of the previous example. Instead of having two species of tiles or supertiles, we have 
three, with the fusion rules 

P n (a) = (P n _ 1 (a)) 10 "P„_ 1 (6)P n _ 1 (c)(P n _ 1 (a)) 10 ", 
Pn(b) = (P t - 1 (6)) 10 "P n _ 1 (a)P n _ 1 (c)(P n _ 1 (6)) 10 ' 1 , 
P n (c) = (P n _ 1 (a)) 10 ' l P n _ 1 (c)P n _ 1 (c)(P n _ 1 (6)) lon , 

/ 2x 10 n 1 10 n \ 

with transition matrix M n _i 5 „ = ( i 2xio n 10™ ). As before, the measure \i a coming from 
the sequence a, a, . . . is ergodic and describes the patterns in a high-order a supertile, which 
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is rich in a tiles, while the measure //& describes the patterns in a high-order h supertile, which 
is similarly rich in b tiles. The measure fi c from c, c, . . . describes a high-order c supertile, 
which is (essentially) half high-order a supertiles and half high-order b supertiles. In other 
words, /i c = (/i a + /ib)/2 is not ergodic. 

3.5. Unique ergodicity. A system is uniquely ergodic if it has exactly one ergodic probabil- 
ity measure, in which case this measure is the only invariant probability measure whatsoever. 
Tiling dynamical systems are uniquely ergodic when there are uniform patch frequencies that 
can be computed regardless of the tiling (see e.g. Theorem 3.3 of [57]). 

For each n, we say that the family of matrices D n ^ is asymptotically rank 1 if there is a 
vector d n e W n such that the columns of D n ^ all approach d n as iV — > oo. Put another 
way, D U) n is asymptotically rank one if A n consists of a single point. 

Theorem 3.9. // a primitive fusion rule 1Z is van Hove and recognizable, then D n ^ is 
asymptotically rank 1 for every n if and only if the tiling dynamical system (X^,]R d ) is 
uniquely ergodic. 

Proof. By Corollary 3.5, having a unique measure is the same as Aoo being a single point, 
which is equivalent to each A n being a single point. □ 

Corollary 3.10. The tiling dynamical system of a transition-regular fusion rule that is 
recognizable, van Hove and primitive is uniquely ergodic. 

What remains is to find checkable conditions on the transition matrices M n ^ that imply 
that the direction matrices D n ^ are asymptotically rank one. For the n-th transition matrix 
M n _i jTl , let 

, c \ * ( mim,M n _ lin (i,/c) 

(3-5) n = mm — — — 

k \max i M ri _i in (z,fc)_ 

This measures the extent to which the columns of M n _i jn are unbalanced. 

Theorem 3.11. If S n diverges, then 1Z is primitive and for each n the family D n ^ is 

n 

asymptotically rank 1. 

Proof. First we show that the diameter of A nj 7v+i is bounded by (1— 5n+i) times the diameter 
of A n) Ar. Let v Uj n be the sum of the columns of M ni Ar, and let v n> N be the direction of v n> N. 
Let mjv + i 5 j be the smallest entry of the ith column of Mn^n+i (which may be zero). The zth 

column of M^^+i is then mjv+i,i ( : J , plus additional terms, so the ith column of M n ^ + i 

is TnN+\,i v n,N-, phis an additional linear combination of columns of M n ^. This means that 
the direction of the ith column of M n ^ + i is a weighted average of v n> ^ and an unknown 
element of A ni N, with v n N having weight at least 5n+i- Thus the direction of each column 
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of M n ^ N+ i, and hence A UtN+1 lies in the convex set 5jv+i"£n.,iv + (1 ~ ^N+i)^n,N, a set whose 
diameter is (1 — <5at+i) times the diameter of A UjN . 

If S n diverges, then S n is nonzero infinitely often, so the fusion rule is primitive. Fur- 

oo 

thermore, the infinite product J~J (1 — 5k) equals zero. Thus A n has diameter zero, and is 

k=n+l 

a single point. □ 

Corollary 3.12. If 1Z is a strongly primitive, van Hove and recognizable fusion rule whose 
transition matrices M„_i jn have uniformly bounded elements, then (X/^IR^) is uniquely er- 
godic. 

Proof. If the smallest matrix element of M n _i jn is at least 1 and the largest is at most K, 
then each 5 n > 1/K. Thus ^^S n diverges and every D n N is asymptotically rank 1. By 

n 

Theorem 3.9, (Xfc,M d ) is uniquely ergodic. □ 

3.6. Transversals, towers, and rank. In tiling theory, especially the aperiodic order and 
quasicrystal branches, the concept of the transversal is an essential component to many 
arguments. For instance, it is used for computing the C*-algebras and i^-theory as in [36] 
and references therein, and it is used for gap-labelling in [ ] . In ergodic theory, the concept 
of towers and especially the Rohlin Lemma (also called the Kakutani-Rokhlin or Halmos- 
Rokhlin Lemma) is a tool that has been used to great effect (see for instance [46, 44]) . One 
notable result [5] that uses towers and the lemma is that any aperiodic measure-preserving 
transformation on a standard Lebesgue space can be realized as a cutting and stacking 
transformation. Towers are used to define the notion of rank, which is intimately related to 
spectral theory. 

For convenience, we will assume that G = G t and that our supertile sets V n are described 
as follows. We position each prototile in Vq so that the origin is in its interior, and the 
place where the origin sits is called the control point of the prototile. In a prototile that 
has been translated by some element v G M. d we call v the control point of the new tile. 
Each element of V\ is positioned such that the origin is on the control point of one of the 
prototiles it contains, and this point will be considered to be the control point of the 1- 
supertile. Likewise, we situate the elements of V2,Vs, etc. in such a way that the origin 
forms the control point of the n-supertile and lies atop the control point of the T^-i-tile at 
the origin. 

Definition 3.13. The transversal of X-ji is the set of all tilings positioned with the origin 
at the control point of the tile that contains it. If 7Z is recognizable, the n-transversal of 
is the set of all tilings positioned so that the origin is at the control point of the n-supertile 
containing it. 
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(If G r is nontrivial, the situation is only slightly more complicated. For each prototile, 
we fix a preferred orientation. The transversal of a tiling space is the set of tilings with the 
origin at a control point, and with the tile containing the origin in the chosen orientation. 
The n-transversals are defined similarly.) 

The transversal of X-ji has a natural partition into jo disjoint sets, one for each type of 
tile. Each of these can be decomposed into pieces, one for each way that the tile containing 
the origin can sit in a 1-supertile, and this partitioning process can be continued indefinitely. 
The n-transversal of X-ji can be thought of as the transversal of X^, and can likewise be 
partitioned. When the fusion rule has finite local complexity, the transversal is a totally 
disconnected set. The n-transversals will form the base for the nth tower representation. 

In one-dimensional discrete dynamical systems, phase space can be visualized as a stack of 
Borel sets placed one above the other with the transformation taking each set to the one di- 
rectly above it, except the top one, on which the action is not visualized. This representation 
of the system is known as a Rohlin tower. When the action is continuous, multidimensional, 
or by an unusual group, the "towers" no longer resemble physical towers, but the term still 
applies. The concept of Rohlin towers for groups other than Z, and in particular for M. d , is 
investigated in [44, 35, 18, 51] and our definitions are drawn from these. Let (X, 23, u) be 
a probability space acted on by some amenable group G to produce an ergodic dynamical 
system. 

Definition 3.14. (1) Let B C 03 and let F C M. d be relatively compact, and suppose that 
g{B) H h(B) = for any g ^ h G F . In this case we call (B, F) a Rohlin tower with 
base B, shape F, and levels g(B). 

(2) A tower system is a finite list of towers T = (Bi, Fx), (B n , F n ) such that all levels 
are pairwise disjoint. 

(3) The support of a tower system is the union of its levels and the residual set is the 
complement of the support in X . 

(4) A sequence T k of tower systems is said to converge to 03 if for every Borel set A G 53 
and every e > there is an N such that for all k > N there is a union of levels of 
T k whose symmetric difference from A measures less than e. 

Recognizable fusion tiling dynamical systems come automatically equipped with tower 
systems that converge to the Borel a-algebra 03 (X^.). The nth tower system will have one 
tower for each prototile P n (j) G V n , for a total of j n towers. The base of the jth tower is 
the set B n (j) of all tilings in the n-transversal that have a copy of P„,(j) with its control 
point at the origin. The shape of the jih tower, denoted F n (j), depends on whether G r is 
trivial. If so, then the shape is the set of all v G G t such that the origin is in the interior of 
P n (j) + v. This shape is, up to sign, the same as the interior of P n (j) itself, thus earning 
the name "shape". If instead G r is nontrivial, the tower construction must be modified to 
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accomodate tilings that have discrete rotational symmetry. In this case F n (j) is the set of 
all g G G for which the origin is in g(P n (j)), provided P n (j) has no symmetry. If it does, we 
must restrict the rotational portion of F n (J) to keep the levels disjoint. 

In ergodic theory an important idea is that of rank. A dynamical system is said to have 
rank r if for every e > there is a tower system (P>i, Fi), (B r , F r ) that approximates all 
elements of 23 up to measure e, where r is the smallest integer for which this is possible. It 
is well-known for substitution sequences and self-affine tilings that the rank is bounded by 
the size of the alphabet or the number of prototiles, since every tile type gets its own tower 
for each application of the substitution. In the case of fusion tilings, the situation is only 
slightly more complicated and we can say 

(3.6) rank(X7£,G) < liminfjn 

n— >oo 

In general, rank bounds spectral multiplicity. 

3.7. Groups other than M. d . For much of this section we have assumed that G = M. d . 
However, the results can readily be adapted to tiling spaces that involve other groups. In 
this section we indicate what changes have to be made when G t is a proper subgroup of M. d , 
when G r is nontrivial, or both. 

In general, G t is the product of two groups, namely a continuous translation in a subspace 
E of M. d , and a discrete lattice L in the orthogonal complement of E. In place of Lebesgue 
measure on R rf , the measure on Gt is the product of Lebesgue measure on E and counting 
measure on L. Frequencies are defined as before as occurrences per unit volume in Gt- In 
fact, the ergodic theorem and Rohlin towers were first developed for discrete group actions 
and only later extended to continuous groups. 

Having G r nontrivial is more of a complication, especially if G r is continuous, as with the 
pinwheel tiling. The ergodic theorem still applies, since we can first average over G r and 
then average over G t , but the G-orbit of a tiling is can no longer be identified with Euclidean 
M. d . When G r is continuous, the frequency of a patch is no longer "number per unit volume", 
but is "number per unit volume per unit angle", and may depend on angle. 

If the group G' that defines our dynamics is the same as the group G used to construct the 
tiling, then invariant measures are parametrized exactly as before, by Aoo, or equivalently by 
sequences of well-defined supertile frequencies. The only difference is that p n (i) is the sum 
or integral over angle of the frequency of the supertile P n (i). That is, it counts the average 
number per unit area of P n (i)'s appearing in any orientation. 

If G' is different from G, then we must distinguish between the G-invariant measures, which 
are parametrized by A^, and the G'-invariant measures, which may not be. Determining 
whether every G'-invariant measure is G r -invariant is a separate computation. 
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4. Spectral theory, entropy, and mixing 

A vector a e R d is a topological eigenvalue of translation if there is a continuous map 
/ : Xft — > S , where S 1 is the unit circle in C, such that, for every T £ Xn and every 



The map / is called a topological eigenfunction. Measurable eigenvalues and eigenfunctions 
are defined similarly, only for / measurable rather than continuous. Of course, since con- 
tinuous functions are measurable, every topological eigenvalue is a measurable eigenvalue. 
Given a translation-invariant measure /x on X-n, one can also ask how translations act on 
L 2 (Xn, /jl). The pure point part of the spectral measure of the translation operators is closely 
related to the set of measurable eigenvalues. is said to have "pure point spectrum" 

if the span of the eigenfunctions is dense in L 2 (Xn,ix). 

The set of topological eigenvalues, the set of measurable eigenvalues, and the spectral 
measure of the translation operators are all G>-invariant. If g £ G r and / is an eigenfunction 
with eigenvalue a, then / og is an eigenfunction with eigenvalue g _1 (df). If G r is continuous, 
then this means that the number of nontrivial eigenvalues is either uncountable or zero. 
The first is impossible, as L 2 (Xn,fi) is separable, and as every topological eigenvalue is 
a measurable eigenvalue. The search for (topological or measurable) eigenvalues only has 
meaning, then, when G r is discrete, in which case we can increase our prototile set by 
counting each orientation separately and take G r trivial. 

Standing assumption (Chapter 4): For the remainder of this section, we assume 



A measurable dynamical system is said to be weakly mixing if there are no nontrivial mea- 
surable eigenvalues. A topological dynamical system is topologically weakly mixing if there 
are no nontrivial topological eigenvalues. For primitive substitution tiling spaces there is no 
distinction between the two sorts of weak mixing, as it has been proven [33, 15] that every 
measurable eigenfunction (with respect to the unique invariant measure) can be represented 
by a continuous function. For fusion tilings, the situation is more subtle. In Theorem 4.1 we 
develop necessary and sufficient conditions for a vector to be a topological eigenvalue. This 
theorem is similar to Theorem 3.1 of [13], and of earlier 1-dimensional results of [ ]. The 
key differences are that we work with M. d actions rather than Z d actions and that we do not 
assume linear repetitivity. 

Unlike the substitution situation it is possible for a fusion tiling space to have a measurable 
eigenvalue that is not a topological eigenvalue. In example 4.4 we exhibit a fusion tiling space 
that has pure point measurable spectrum but that is topologically weakly mixing. 

4.1. Topological eigenvalues. For self-affine substitution tilings there are well-established 
criteria for a vector being an (topological or measurable) eigenvalue [57]. Given a substitution 






G = R d . 
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with stretching map L, there is a finite list of vectors Vi such that a is an eigenvalue if and 
only if 

(4.2) a • L n (#;)->• (modi) 

for each i. 

Our first task is to construct an analogous criterion for topological eigenvalues of fusion 
tilings. Assuming strong primitivity, each (n + 2)-supertile contains multiple copies of each 
n-supertile (at least one per (n + 1) supertile). Let V n be the set of relative positions of two 
n-supertiles, of the same type, within an (n + 2)-supertile. This is a finite set, since there 
are only finitely many kinds of (n + 2)-supertiles and each (n + 2)-supertile contains only 
finitely many n-supertiles. We call the elements of V n return vectors. For each a G M. d , let 
T] n ((x) — max I exp(27rza • v) — 11. 

Theorem 4.1. Let 71 be a strongly primitive and recognizable van Hove fusion rule with 
finite local complexity. A vector a G M. d is a topological eigenvalue if and only if 2^ Vn{&) 

n 

converges. 

For primitive substitution tilings, V n+1 = LV n and T] n (a) either goes to zero exponentially 
fast, or does not go to zero at all [ ]. In such cases, the convergence of f] n (a) is equivalent 

n 

to r) n (a) — > 0, which is equivalent to the criterion (4.2), where the vectors range over V°. 

Proof. Since X-n is minimal, a continuous eigenfunction with a given eigenvalue is determined 
by its value on a single tiling T. Fix T G Xn and a G M. d and define /(T) = I. For each 
x G IR d let /(T — x) — exp(27rza ■ x). If this function is continuous on the orbit of T, then 
it extends to an eigenfunction on all of X^. If it is not continuous, then a cannot be a 
topological eigenvalue. 

Suppose that ^^r/ n (df) converges. We will show that / is continuous on the orbit of T. 

oo 

Choose e > and pick n large enough that ^^7/ n (a) < e. We will show that if T — x and 

k=n 

T — y agree to the nth recognizability radius p n , then /(T — x) and /(T — y) are within e. 
The following lemma states that y — x can be expressed as a sum of return vectors. 

Lemma 4.2. Suppose that x and y are corresponding points in n-supertiles of the same type 

N-2 

within the same N -supertile, with N > n + 2. Then y — x can be written as Vk, where 

k=n 

Vk G V k . 

Proof of lemma. For each n, we work by induction on N. The base case N = n + 2 follows 
from the definition of V n . Now suppose the lemma is true for N = No, and suppose that 
x and y sit in the same (No + l)-supertile. The point x sits in an (No — 1) supertile S x , 
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say of type i, and y sits in an iVo-supertile, say of type j. By strong primitivity, there is an 
(iVo — 1) supertile S y of type i in the iV -supertile that contains y. Let zbe the point in S y 
corresponding to where x sits in S x . (See Figure 1.) Then z — x is a return vector from S x 




(N +l) '-supertile 



Figure 1. The induction step identifies a return between the n-supertiles 
(shown as shaded triangles) inside their (N + l)-supertile. 



to S y and so equals £>V -i £ V^ 1 . Meanwhile, y and z sit in the same kind of n-supertile 

N -2 JVo-1 



as 



within the same A^-supertile, so y — z = v^. This means that y — x = v^, 

k=n k=l 

desired. □ 

N-2 

If x and y lie in the same N-th order supertile, the lemma implies that y — x = Vk, so 

k=n 

\f(T-y)-f(T-x)\ = \exp(2ma-(y-x))-l\ 

N-2 N-2 
< | exp(27rzct • V k ) — 1| < T]n(&) < e - 

k=n k=n 

Even if a; and y do not lie in the same iV-supertile of T for any N, it is still true that any 
patch containing x and y is congruent to a patch that lies within an iV-supertile, so y — x 

N-2 

still takes the form Vk for some N and we still obtain that |/(T — y) — /(T — x)\ < e. 

fc=n 

This estimate proves that / is continuous on the orbit of T. By minimality, it extends to a 
continuous eigenfunction on all of X-r. This proves half of Theorem 4.1. 

For the converse, suppose that ^^7y n (cJ) diverges. Then there exists a subsequence 

n 

r] n+3 fc(a) that also diverges. We have the following lemma, that states that any finite sum 

k 

of return vectors separated by three levels appears in X-ji as the return of some n-supertile 
to itself. 
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Lemma 4.3. Pick N such that N + 1 — n is divisible by 3, and for k = n, n + 3, . . . ,N — 2 
pick 4 6 V fc . Let v = v n + v n+3 + ■ ■ ■ + vn-2- For every such set of choices, there exists an 
N-supertile containing two n-supertiles of the same type, such that the relative position of 
the two n-supertiles is v. 

Proof. Again we work by induction on iV. If iV = 2, then this follows from the definition of 
V n . Now suppose it is true for N = N Q , and we shall attempt to prove it for N = N + 3. By 
the inductive hypothesis, there exist points x and y in corresponding n-supertiles within 
the same iV -supertile Si such that y — x = v = v n + v n+3 + • • • + Vn -2, and there exist 
two (A^o + l)-supertiles 5*2 and S3, of the same type and with relative position vn +i, within 
an (iVo + 3) supertile. (See Figure 2). By primitivity, S2 and S3 each contain copies of Si (in 




(Nq+ 3) '-supertile 

Figure 2. The induction step gives a return vector y — x between two copies 
of Si, shown shaded inside an (No + 3)-supertile. 

corresponding positions). Let x be the point corresponding to x in the copy of Si inside S 2 , 
and let y be the point corresponding to yo in the copy of Si inside S3. Then y — x = v. □ 

Thus for any e and for infinitely many values of n, one can find vectors v n , v n+3 , • • • , £^-2 
with Vk G V k , such that | exp(27ria • / ]vk) ~M > 2e. By retricting to a subsequence we can 
assume that the complex numbers exp(27ria • Vk) are either all in the first quadrant or all in 
the fourth quadrant, and that | exp(27rza • ^^Vk) — 1| > e - By Lemma 4.3 there then exist, 

for n arbitrarily large, two n-supertiles of the same type with relative position V^u*- Pick 
x and y to be corresponding points of these supertiles in T, such that a big ball around x 
and y lies entirely in the supertile. (This is possible since the supertiles form a van Hove 
sequence.) Then /(T — x) and f(T — y) differ in phase by at least e, so our purported 
eigenfunction is not continuous. □ 

4.2. Measurable eigenvalues. In this section we provide an example that has pure discrete 
spectrum from a measurable standpoint while being weakly mixing from a topological one. 
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Example 4.4. The scrambled Fibonacci tiling. We consider four fusions, denoted by the 
letters J 7 , A, £ and S, with the last being the "scrambled Fibonacci" fusion whose tiling 
space has the desired properties. All use the prototile set {a, b} where the length of a is the 
golden mean and the length of b is 1. 

The first fusion rule is the usual Fibonacci rule J 7 , which is prototile- and transition-regular 
with (n + l)-supertiles given by F n+ i(a) = F n (a)F n (b), F n+ i(b) = F n (a). This fusion rule 
generates the self-similar Fibonacci tiling space Xjr, which is known to have measurable and 

topological pure point spectrum with eigenvalue set —=Z[(p]. Importantly for our calcula- 

V5 

tions, the Euclidean length of F n _i(a) and F n (b) is 0™, which deviates from an integer by 
±0~ n , and differs from a multiple of 1/0 by ±(f)~( n+2 \ The transition matrix is M = ( \ \). 
To make the second fusion rule, called "accelerated" Fibonacci, we first fix some increasing 
sequence of positive integers {iV(n)}^L 1 , where we assume N(n) — iV(n — 1) > 2. We define A 
to be the induced fusion rule on N(n) levels so that A n (a) = i*jv(n)( a ) an d A n (b) = F^^(b). 
The lengths of the a and b n-supertiles for A are Ar( ™) +1 and (j) N ^ n \ respectively. This fusion 
rule is prototile-regular but not transition-regular, since now the transition matrices M n are 
given by M w v ' . 

The third fusion, which we call "exceptional", is derived from the accelerated rule by 
introducing a third supertile type on all odd levels and using it to introduce a relatively 
small defect in the next (even) level. On both the odd and even levels of £, the supertiles 
E n (a) and E n (b) are constructed with the same populations of prototiles and (n— 1) supertiles 
as A n (a) and A n (b). When n is odd, and only when n is odd, there is an additional tile 
E n (e) with the same population of (n — l)-supertiles as E n {b). The fusion rules for £ n are 
as follows: 

When n is odd, E n (a) and E n {b) are built from £ n ~i in exactly the same way that the 
^4. n -supertiles are built from A n -i- The exceptional supertile E n (e) is obtained from E n {b) 
by permuting the (n — l)-supertiles so that all of the _E n _i(a)-supertiles come before any of 
the i? n _i(6)-supertiles. When n is even, E n (a) and E n {b) are built from £ n -\ in exactly the 
same way that the ^4 n -supertiles are built from A n -i, except that one copy of E n _i(b) is 
replaced with E n _i(e). 

We can make the fusion rule prototile-regular by taking the induced fusion of £ on even 
levels. We call this last fusion rule the scrambled Fibonacci S, but most of our proofs 
center on the equivalent space Xg. By controlling the sequence N(n) we can change spectral 
properties of the scrambled Fibonacci fusion. 

The space Xj is well known to be recognizable, and the recognizability of Xjs, is similar. 
The same patterns that allow us to recognize supertiles in Xj also work (with small modi- 
fications) in X £ and X s . We can thus freely speak of the (unique) n-supertile containing a 
particular tile. 
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Proposition 4.5. If N{2n) - N{2n - 1) goes to infinity fast enough that ^ c^ 2 ™" 1 )-^ 2 ™) 

n 

converges, then all four fusion spaces have pure point measurable spectrum with eigenvalues 
i=(Z + 0Z) = -^Z[0]. 

Proof. We will show that the four spaces Xj?, X_4, Xg and X$ are all measurably conjugate. 

Then, since Xj is well known to have pure point spectrum with eigenvalue set — pZ[</>], the 

V5 

others must as well. Since Xjr and X4 are manifestly the same, and since Xg — X§, we 
need only show that X4 and Xs are measurably conjugate. 

In the tilings in Xg, we call a supertile of any level exceptional if it lies in an E m (e)- 
supertile of some level. Note that E 2m+2 (a) and E 2m+2 {b) each contain only one E 2m+ i(e) 
supertile and of order (j)N^m+2)-N(2m+i) supertiles of type E 2m+1 (a) and E 2m+1 (b). The 

fraction of (2n)-supertiles that are exceptional in X| n is thus bounded by a constant times 

00 

e n = 0-(JV(2m+2)-iV(2m+i))^ which by assumption goes to zero as n — > 00. 

m=n 

Suppose T is a tiling in Xg. If the origin lies in an unexceptional supertile of some level 
n, and hence also at levels n + 1, n + 2, etc., and if the union of these supertiles is the 
entire line, 3 then we can convert this to a tiling in X^ by replacing each unexceptional 
^-supertile containing the origin with the corresponding ^4-supertile. From the definitions 
of the supertiles, this operation 01111 + 1-supertiles is consistent with the operation on n- 
supertiles. 

The measure of the tilings for which the origin is in an exceptional n-supertile is bounded 
by a constant times e n , and so goes to zero as n — > 00. Thus, with probability 1, the origin 
lies in an unexceptional supertile of some level. Likewise, with probability 1, the union of 
the supertiles containing the origin is all of 1R. Thus we have a map from Xg to X4 that is 
defined except on a set of measure zero. This map is readily seen to preserve measure and 
to commute with translation. □ 

Proposition 4.6. // lim N{2n+1) — 2N(2n) = +00, then Xs is topologically weakly mixing. 

n—>oo 

Proof. First we show that elements of ^pZ[0] cannot be topological eigenvalues. Then we 

V5 

show that real numbers that are not of this form cannot be topological eigenvalues. 

The supertiles S n (a) = E 2n (a) and S n (b) = E 2n {b) have length 0^ 2n ) +1 and (f) N{2n \ 

respectively. For a e — pZ[</>], | exp(27ria|S' n (a)|) — 1| and | exp(27rza;|5' n (6)|) — 1| are bounded 

V5 

above and below by constants (depending on a) times cf)~ N<y2n \ Each supertile E 2n+ i(e) 
contains the Fibonacci number fN(2n+i)-N{2n) consecutive copies of E 2n (b), since that is how 
many iV(2n)-supertiles of type b there are in the N{2n + l)-supertile i 7 jv(2n+i)(&)- We thus 

3 Note that this condition is translation-invariant, as every point in T would then lie in a sequence of 
unexceptional supertiles whose union is the entire line. 
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find at least that many consecutive copies of E 2n (b) in S n+2 ( a ) and S n+ 2(b), so there exist 
vectors Vk = k\S n (b)\ between n-supertiles of the same type in the same (n + 2)-supertile, 
where k is any positive integer up to fN(2n+i)-N(2n)- Since | exp(27rza|5' n (6)|) — 1| is bounded 
below by a constant times <p~ N ^ 2n \ and since fN(2n+i)-N(2n) is of order ( j ) N ( 2n + 1 )- N ( 2n ) ? anc [ 
since 7V ( 2n+1 )- 2iV ( 2 ") grows without bound, there are k for which for which exp(2niav k) is 
not close to 1. In fact, by taking n sufficiently large and picking k appropriately, we can get 
exp(27riat>fc) to be as close as we want to any number on the unit circle. 

On the other hand, if a is not in —^=Z [</>], then by Pisot's theorem, exp(27ri0 n a) does not 

V5 

approach 1 as n — > oo. Since for arbitrarily large patches P there exist return vectors of 
length <fr n for n sufficiently large, a cannot be a topological eigenvalue. □ 

It is simple to construct sequences N(n) that meet the conditions of both Propositions 
(4.5) and (4.6). For instance, we could take N(n) = 3 n . Thus there exist fusion tilings that 
are topologically weakly mixing but are measurably pure point. 



4.3. Pure point spectrum. An important and widely studied problem in substitution se- 
quences and substitution tilings is determining when the (measure-theoretic) tiling dynamics 
have pure point spectrum. A key tool is Dekking's coincidence criterion [19], first developed 
for 1-dimensional substitutions of constant length and later extended to arbitrary substitu- 
tions, with generalizations in higher dimensions such as Solomyak's overlap algorithm [57]. 
In this section we explore the extent to which the analog of Dekking's criterion determines 
spectral type for fusions. The differences between substitutions and fusions are already 
apparent in the simplest category, namely one dimensional fusions of constant length. 

We say a 1-dimensional fusion (or substitution) has constant length if, for each n, all of 
the n-supertiles P n (i) have the same size. This implies that tiles all have the same length 
and that, for fixed n, each n-supertile contains the same number L n of (n — l)-supertiles. 
The fusion is coincident if, for each n, there exists an N such that any two A-supertiles 
agree on at least one n-supertile. The fusion is coincident with finite waiting if there exists a 
fixed integer k such that N = n + k works for every n. For substitution tilings, coincidence 
is equivalent to coincidence with finite waiting, but for fusions it is not. 

To each fusion of constant length we associate a solenoid S-r, obtained as the inverse 
limit of the circle IR/Z under a series of maps, with the n-th map being multiplication by 
L n . S-r is a topological factor of X-r, with a point in S-r describing where the origin lies 
in a tile, a 1-supertile, a 2-supertile, etc, but not generally determining which type of n- 
supertile the origin sits in. There is a natural translational action on S-r, and the span of 
the eigenfunctions of this action is dense in L 2 (S-r). If the factor map from X-r to S-r is 
a measurable conjugacy, or equivalently if there is a set of full measure on X-r where the 
factor map is 1:1, then X-r has pure point spectrum. If the factor map is not a conjugacy, 
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and if every eigenfunction on Xr is obtained from an eigenfunction on Sr, 4 then Xr does 
not have pure point spectrum. 

For substitutions of constant length, the situation is clear-cut: 

Theorem 4.7 ([19]). A one dimensional tiling space obtained from a primitive and recog- 
nizable substitution of constant length and height one has pure point measurable spectrum if 
and only if it is coincident. 

There are two reasons why this theorem does not apply to general fusions. First, a 
coincident fusion may not be uniquely ergodic. For each ergodic measure, the question isn't 
whether a generic point in Sr corresponds to a single tiling, but whether it corresponds 
to a single tiling in a suitably chosen set of full measure. Second, a coincidence, or even a 
coincidence with finite waiting, only implies that supertiles agree somewhere. Unless we have 
some control over the transition matrices, we cannot conclude that high-order supertiles agree 
on a fraction approaching 1 of their length, which is what is needed to obtain a measurable 
conjugacy between Xr and Sr. 

In Example 3.7, the fusion is not coincident, as P n (a) and P n (b) disagree at every site. 
The map from Xr to S-r is nowhere 1:1, being 4:1 over the orbit where there exist two 
infinite-order supertiles, and 2:1 over all other orbits. 

However, for each ergodic measure, Xr does have pure point spectrum. For instance, for 
the ergodic measure that comes from the supertile sequence {a, a, . . .}, the measure of the 
tilings where the origin sits in a supertile P n (a) is exponentially close to 1. With probability 
1, for all sufficiently large n the n-th order supertile containing the origin is of type a. Also 
with probability 1, the infinite-order supertile containing the origin covers the entire line. 
The set of tilings with both these properties has full measure, and the factor map to Sr is 
1:1 on this set. 

Example 4.8. To see how coincidence with finite waiting is insufficient to prove pure point 
spectrum we make a fusion based on the substitution cr(6) = bc 5 b 4 , cr(c) = cb 5 c 4 . Repeated 
substitution produces words c n {b) and o""(c); by abusing notation we write cr n (P n _i(fe)) and 
a n (P n -i(c)) to mean the fusion of {n — l)-supertiles of types b and c in the order given by 
the letters of a n {b) and o~ n (c). We introduce coincidence with finite waiting by defining the 
fusion rule to be 

P n (b) = P n _ 1 (6) ( 7 n (P n _ 1 (fo))P n _ 1 (c), P n (c) = P n ^(b)a n (P n ^(c))P n ^(c). 
The transition matrix M n _i 5 „ = (| |) n + ({ {) = (5 x lO™" 1 + 1) ({ {) has rank 1, so 

n 

the system is uniquely ergodic. The length of an n-supertile is ]^J(10 J + 2), and P n (b) and 

i=i 

4 This is connected to the height of a substitution or fusion. If a substitution has height one, then all 
eigenvalues of X-r are eigenvalues of S-r [45]. One can similarly define a notion of height for fusions. 
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P n (c) differ on lO- 7 (n — l)-supertiles, implying that they differ on J^J 10 J tiles. Thus P n (b) 

and P n (c) agree on a positive fraction of their tiles, namely 1 — J| — — -. As n — > oo, 

3=1 3 + 

this fraction increases but does not approach 1. This implies that one can find disjoint 
measurable sets of positive measure that map to the same set on the solenoid. Any function 
that distinguishes between these sets cannot be approximated by a function on the solenoid, 
so the span of the eigenfunctions is not dense and the spectrum is not pure point. For an 
example of such a function, let /(T) equal 1 if the origin is in a b tile and if the origin is 
in a c tile. 

Example 4.9. It is also possible for different ergodic measures for the same fusion to have 
different spectral types. Consider the 1-dimensional non-primitive substitution 

a(a) = a 10 , a(b) = bc 5 b\ a(c) = cb 5 c 4 

of constant length 10. Next consider a 1-dimensional fusion tilings with three prototiles 
a, b, c, each of unit length. Using the same abuse of notation as in the previous example we 
define the fusion rule as 

P n (a) = P n _ 1 (a) ( 7"(P n _ 1 (a))P n _ 1 (6)P n _ 1 (c) 
P n (b) = P n _ 1 (a) < r n (P n _i(6))P n _i(6)P n _i(c), 
Pn(c) = P n _ 1 (a)a n (P n _ 1 (c))P n _ 1 (6)P n _ 1 (c), 

This fusion is coincident with waiting time 1, in that all n-supertiles begin with P n _i(a) and 
end with P n _ 1 (fe)P n _ 1 (c). However, that is only 3 out of L n = 10™ + 3 (n — l)-supertiles, 
and the n-supertiles disagree on the rest! For large n, the supertiles P n (a), P n (b) and P n (c) 
disagree at roughly 70% of their tiles, 97% of their 1-supertiles, 99.7% of their 2-supertiles, 
etc. 

The transition matrices 

/I0 n + 1 1 1 \ 

M n _ 1>n = 1 5 • 10"" 1 + 15- 10"- 1 + 1 
\ 1 5-10 n_1 + l 5-10 n_1 + l/ 

have rank 2. There are two ergodic measures, one coming from the supertile sequence 
{a, a, . . .} and the other coming from an arbitrary sequence of fe's and c's. 

When we take the ergodic measure from the sequence {a, a, . . .}, X-ji is measurably con- 
jugate to the solenoid S-r and has pure point spectrum. When we take the other ergodic 
measure, however, there is a set of full measure where, for all sufficiently large n, the origin 
is either in P n (b) or P n (c), but is not in the two right-most n — 1 supertiles within P n (b or c). 
This set admits a measure-preserving involution where, for all sufficiently large n, the su- 
pertiles P n (b) containing the origin are replaced by P„(c) and vice- versa. On any set of full 
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measure, the factor map is (almost everywhere) at least 2:1, and the tiling dynamical system 
does not have pure point spectrum. 

In other words, almost every point of the solenoid corresponds to three tilings. One set of 
preimages has full measure with respect to the {a, a, . . .} ergodic measure, while the other 
two preimage sets have full measure with respect to the other ergodic measure. Since the 
first ergodic measure only "sees" one preimage, it has pure point spectrum. Since the other 
measure "sees" two preimages, it does not have pure point spectrum. 

To get pure point spectrum from coincidence, we must control the transition matrices. 

Theorem 4.10. Let 71 be a primitive, recognizable, prototile-regular, 1- dimensional fusion of 
constant length. If the fusion is coincident with finite waiting, and if the transition matrices 
M n _i tn are uniformly bounded, then Xr is uniquely ergodic and has pure point spectrum. 

Proof. Suppose that there are J species of prototiles, that the fusion is coincident with 

waiting k, and that M n _ l n (z,j) < C for all n,i,j. Then L n < CJ. Any two n/c-supertiles 

agree on at least one (n — l)fc-supertile, at least one (n — 2)fc-supertile in each remaining 

(n — l)/c-supertile, at least one (n — 3)fc-supertile in each remaining (n — 2)£>supertile, etc. 

C k J k -V v 



This means that any two nfc-supertiles agree on at least a fraction 1 — — —r — r — of their 

\ C k J k J 

tiles, a fraction that approaches 1 as n — > oo. In particular, the density of tiles (and likewise, 
of n-supertiles for any fixed n) is asymptotically the same in all iV-supertiles as iV — > oo, 
implying unique ergodicity. 

A point in S-r thus determines all but a set of density zero of the tiles in the infinite-order 
supertile containing the origin. The probability of there being an undetermined tile in any 
fixed bounded region is thus zero. Since the real line is a countable union of bounded regions, 
and since the probability of having two infinite-order supertiles in a single tiling is also zero, 
almost every point in S-r corresponds to a tiling with no undetermined tiles. Thus the factor 
map Xr S-r is a measurable conjugacy, and Xr has pure point spectrum. □ 

Theorem 4.10, while modest in scope, is typical of the theorems that can be proven about 
fusions that are not of constant length, or that are not 1-dimensional. Given any concidence- 
based test for pure point spectrum in the category of substitution tilings (e.g., the balanced 
pair algorithm or the overlap algorithm), one can construct an analogous test for fusions. 
However, a positive result from such a test will only demonstrate pure point spectrum if one 
can also show that the coincidences happen frequently enough. This requires estimates both 
on how long one must wait for a coincidence, and on how much the system has grown in the 
process. 



4.4. Entropy. Standard results in symbolic and tiling substitution dynamics say that such 
systems cannot have positive entropy. The obstruction is that the transitions from level 
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to level do not contain much 'new' information. This continues to be the case for fusion 
tilings when one assumes that both the number and shapes of supertiles remain fairly well- 
controlled. We give a simple example of a minimal and uniquely ergodic fusion rule with 
positive entropy, and then we prove a sufficient condition for a fusion space to have zero 
entropy. 

Configurational entropy is based on counting configurations, and for this we need G to be 
discrete. We therefore assume that G = Z d , so we are essentially dealing with subshifts. Let 
be the number of configurations that can appear in a d- dimensional cube of side n (this 
is the complexity function). The configurational entropy is 

log #n 



lim 



n d 



For subshifts, configurational entropy is known to be the same as topological entropy. 

Positive entropy implies that there is a lot of randomness in the system, while unique 
ergodicity means that all patterns appear with well-defined frequencies. These ideas might 
seem to be in conflict, but Jewett [34] and Krieger [■ J] showed that uniquely ergodic dynami- 
cal systems can exhibit a very wide range of dynamical properties, and in particular can have 
positive topological entropy. The following example is in the spirit of their construction. 

Example 4.11. A strictly ergodic fusion rule with positive entropy. We construct a fusion 
rule 1Z with P = {a, b} recursively. Let V\ be all words of length 3 in which each letter 
appears at least 1 time but no more than 2 times; we have j± = 6 distinct 1-supertiles. Now 

let Vi be all fusions of = 54 1-supertiles in which each supertile appears between ji and 

2ji times. The expected number of 1-supertiles in any fusion of of them is so we 
are including the highest-probability fusions in our set Vi- 

3j2 

In general, let j n be the number of n-supertiles and let V n +i be all fusions of n- 
supertiles in which each n-supertile appears between j n and 2j n times. Since having more 

3j2 

than 2j n or fewer than j n occurrences in a span of size — - is already highly improbable, 
restricting to these (n + l)-supertiles only reduces the number of configurations slightly, and 
the system so constructed has positive entropy. The transition matrices are enormous and 
grow super-exponentially but always have all nonzero entries, making the system strongly 
primitive and hence minimal. Moreover, the constant 5 n used in equation (3.5) (to measure 

Jn 

how balanced the columns of the transition matrices are) is always = 1/2, so the tiling 
space is uniquely ergodic. 

The fusion rule 7Z is not recognizable, but we can build a recognizable fusion rule 1Z' from 
1Z as in Example 2.3. Since the entropy of the factor X n is bounded by the entropy of X-j^j, 
Xn> has positive entropy. It is easy to check that the addition of subscripts does not affect 
unique ergodicity. 
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This example involved the number j n of n-supertile types growing exponentially with the 
size of the supertiles. If the growth is slower than exponential, and if the shapes of the 
supertiles are not too distorted, then the system will have zero entropy. 

Proposition 4.12. Let d n be the diameter of the largest n-supertile, let j n be the number of 
n-supertile types and suppose that there exist constants (3, K such that each cube of side (3d n 

touches at most K n-supertiles. If lim — = 0, then the configurational entropy of X-r, 

ra->oo d^ 

is zero. 

Proof. To determine the configuration in a cube of side f3d n) one must specify the kinds 
of n-supertiles that intersect that cube, and also specify the locations of those supertiles. 
There are at most choices for the first, and at most V choices for the second, where V 
is the volume of the largest n-supertile, which is bounded by d n . Thus the log of the number 

of configurations, divided by the volume of the cube, is bounded by S_n 

p d d d n 

which goes to zero as n — > oo. □ 

The upshot of Proposition 4.12 is that positive entropy either requires the number of 
supertiles to grow exponentially with volume, or for the shapes and relative sizes of supertiles 
to be so wild, and for the ways that supertiles fit together be so varied, that there are many 
ways for supertiles to fit together. The geometric issues do not apply to dimension 1, where 
supertiles simply concatenate, but could in principle apply in dimensions 2 or more. However, 
we know of no examples where positive entropy is achieved without j n growing exponentially 
fast. 



4.5. Strong mixing. A (measurable) dynamical system is strongly mixing if for any pair 
of measurable sets A, B, and for any sequence of vectors v n tending to infinity, lim fi(A fl 
(B — v n )) = ji{A)ji{B). The dynamical systems of primitive substitution sequences and self- 
similar tilings are never strongly mixing [20, 57]. Because of the rigidity of the substitution 
process, knowing the location of one copy of a patch gives a higher probability that it will be 
seen again at certain intervals. However, there are "staircase" cut-and-stack transformations 
in one and several dimensions that have been shown to be strongly mixing [1, 2], thus it is 
possible to have strongly mixing fusion tiling systems. As in the case of entropy, this is only 
possible when the system has increasing complexity at higher levels of the hierarchy. 

In this section we establish sufficient criteria for fusion tilings not to be strongly mixing. 
These criteria involve both uniform bounds on the number of supertiles and on the transition 
matrices, and are not necessary criteria. For instance, the accelerated Fibonacci fusion 
discussed in Example 4.4 does not have bounded matrices, but is essentially the same as 
ordinary Fibonacci and is not strongly mixing. 
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Theorem 4.13. The dynamical system of a strongly primitive van Hove fusion rule with a 
constant number of supertiles at each level and bounded transition matrices, and with group 
G = W 1 , cannot be strongly mixing. 

Proof. Our proof is an adaptation of Solomyak's [57], which in turn is an adaptation of 
Dekking and Keane's [ ]. By Corollary 3.12, X-r is uniquely ergodic, so for any patch 
P, freq a {P) can be computed from the actual frequency of P in any increasing sequence 
of supertiles. We will find a patch P and a sequence of vectors v n , tending to infinity, 
such that the frequency of P U (P + v n ) is bounded away from zero. Then, supposing that 
freq n (P) = 5 and freq u (P U (P + v n )) > e, we pick a set U C M. d whose volume is less than 
— rj, and which is small enough that fi(Xp^u) = freq a (P)Vol(U). Let A — B — Xpjj. Since 
X PU (p+v n ),u cAn(B-t? n ),we have 

(4.3) fi{A n (B - v n )) > ti(X PuP+ ^u) > eVol(U) > 25 2 Vol(U) 2 = 2fji(A)(i{B), 

so fi(A fl B — v n ) cannot approach /j i (A)/j i (B) as n — > oo. 

To find the vectors v, we suppose the number of supertiles at each level is the constant J. 
Since X^ can be expressed as the union of J cylinder sets defined by which n-supertile is at 
the origin, it must be that at least one of those cylinder sets has measure at least 1/J. For 
each n, choose l n G {1,2, ...J} corresponding to an n-supertile with this property, so that 
Vol(P n (l n ))p n (l n ) > 1/J, where p n is the supertile frequency vector. Now choose v n e V™ 
to be a return vector for P n {l n ) inside P n+2 (/ n +2), as in Theorem 4.1. Because of strong 
primitivity and the fact that our transition matrices are uniformly bounded, we can find 

a 5' > for which J n [ 7 n ^\, > 5' for all n. (Specifically, if each n-supertile contains 

Vol(P n+2 {ln+2)) 

at most K (n — l)-supertiles, then the ratio of volume between the largest and smallest 
n-supertile is at most K, and Vol(P n (l n )) /Vol(P n+ 2(l n +2)) > V-^ 3 -) 

The patch P is arbitrary. By choosing n large we can make — — - — tttv~ arbitrarily 

Vol(P n (l n )) 

close to freq^(P), and hence greater then a fixed constant freq fJi (P)/2 for all n. The reader 
can refer to Figure 3 to see that #(P U P + v n in P n +2^n+2)) > #(P i n Pniln))- Since the 
fraction of area from the supertiles P n (l n ) is at least 1/J, this implies that the frequency of 

r, /r, i freqJP)Vol(P n (l n )) , , 5'freqJP) 

P U (P + v n ) is at least J — ^ \' — n " , hence at least J 7^ ; . □ 

2JV 0l{P n+ 2{ln+2)) 2.J 

Proposition 4.14. The dynamical system of a strongly primitive van Hove fusion rule with 
a constant number of supertiles at each level and bounded transition matrices, and with group 
G = 7L d , cannot be strongly mixing. 

Proof. The previous proof does not apply to 7L d actions because we cannot choose U arbi- 
trarily small. However, we have already shown that for any patch P and any sufficiently 
large n, there exist large v with freq^iP U (P + v)) > 5' freq n (P)/2J. We then find a patch 
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Figure 3. Each copy of P in P n {l n ) makes a copy of P U (P + v n ) in P n+2 (/ n +2)- 

P whose frequency is less than 8' /AJ, so that freq^((P U (P + v n )) > 2/reg /i (P) 2 . Taking 
U to consist of one point, this implies that p,(A n (A — v n )) > 2fj,(A) 2 , where A = Xpp. □ 

5. Inverse limit structures, collaring, and cohomology 

In this section we consider topological properties of spaces X n of fusion tilings, including 
their structure as inverse limit spaces, their Cech cohomology groups, and the significance 
of these groups. 

Standing assumptions (Chapter 5): We assume, for simplicity, that G = IR d ; modifi- 
cations for other groups are done exactly as for substitution tilings [43, 55]. We also assume 
that our fusion rules are recognizable and, as always, have finite local complexity. 

Tiling spaces can always be represented as inverse limits of CW complexes [9, 54]. The 
challenge is finding a representation that allows for efficient calculation and for the proving 
of theorems. To this end we present generalizations of the Anderson-Putnam complex [4] 
and of the partial collaring scheme of Barge, Diamond, and their collaborators [6, 7]. 

In all cases, we construct a sequence of spaces and maps 

T J° T J 1 T J 2 T J 3 

1 o i — ip — i 2 — 1 3 — " " - j 

where each approximant Tj describes a region of the tiling, each Tj+i describes a larger region 
of the tiling, and fi : — » Tj is the forgetful map that loses the additional information 
carried in The inverse limit Um(T, /) is the set of infinite sequences (xo,x\, . . .) such 

that each X{ G Tj and each x« = fi(x i+ i). Such a sequence is a set of consistent instructions 
for tiling larger and larger regions of R . If the union of these regions is all of IR d for all 
sequences in the inverse limit, then there is a natural isomorphism between h±m(r, /) and 
Xft. The key is to make sure that every tiling in X-ji can be built up from the approximants 
in a unique way. A common obstruction is when the approximants can build an infinite tiling 
that covers only a portion of M. d . "Border-forcing" fusions, discussed next, do not have this 
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obstruction. Later we will describe the technique of "collaring" to make fusion rules become 
border-forcing. 

5.1. Forcing the border. A fusion rule always tells us how n-supertiles make up the interi- 
ors of larger iV-supertiles. But sometimes the A-supertiles also determine which n-supertiles 
are on their exterior as well. When this happens we say the fusion rule forces the border, 
and we have a natural way to see the space as an inverse limit. 

Definition 5.1. A fusion rule forces the border if for each integer n there exists an N with 
the following property: If S\ and S2 are two N-supertiles of the same type appearing in tilings 
T\ and T2 in X-ji, then the patch of n-supertiles that touch S\ in T\ is equivalent to the patch 
of n-supertiles touching S2 in T2. 

Example 5.2. Compare and contrast: border forcing. The 1-dimensional substitution a — > 
abb, b — > abbb forces the border in that every n + 1-supertile of type a is preceded by an n- 
supertile of type b and followed by an n-supertile of type a, and likewise every n+ 1-supertile 
of type b is also prededed by an n-supertile of type b and followed by an n-supertile of type a. 
By contrast, the substitution a — )■ ab, b — )■ aa does not force the border, since an A-supertile 
of type a can be preceded either by an n supertile of type a or b. 

5.2. The Anderson-Putnam complex. To build r , we start out with one copy of each 
prototile from P - Then, if somewhere in some tiling two prototiles meet, we identify the 
corresponding points on their boundaries. The resulting branched manifold is compact [4]. 
(If we take the periodic tiling of unit squares lined up edge-to-edge, it is easy to see that 
T is the torus.) We build l\ by taking one copy of each supertile from V\ and identifying 
the boundaries whenever they meet as above, and continue making each approximant r„ 
similarly. Put another way, T n for the space X-ji is To for X^. 

There is a natural map from Xn to T n that maps a tiling to the location of the origin 
within its n-supertile. Thus, a point in T n can be viewed as a set of instructions for placing 
an n-supertile containing the origin. Or course, if we know the (n + l)-supertile containing 
the origin, then we necessarily know the n-supertile containing the origin, so the forgetful 
map /„ is well-defined. 

Theorem 5.3. If the recognizable fusion rule 1Z forces the border, then X-ji is homeomorphic 
to the inverse limit h_m(r ra , / n ) of Anderson-Putnam complexes. 

Proof. We will construct a homeomorphism from the inverse limit to X-ji by constructing 
maps from each approximant to partial tilings of lR d . Pick an increasing sequence of integers 
n 1; n 2 , . . . such that all the n-supertiles bounding an nj +1 -supertile are determined by the 
type of the n i+ i-supertile. Our map takes a point xn in to the A-supertile with the 
origin where xn is, together with all of the lower order supertiles that are determined by 
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that iV-supertile. If N > rii, then this includes all the rij^-supertiles touching the N- 
supertile, all the nj_ 2 supertiles touching the rii-i supertiles, all the n;_ 3 supertiles touching 
the rij_2 supertiles, and so on. If x^ is on the boundary of an A-supertile, then there are 
multiple tilings that can come from this process, but they all agree on the nj_i-supertiles 
in all directions around the origin, as well as the lower-order supertiles determined by the 
nj_i-supertiles. In particular, xn determines at least i — 1 layers of supertiles of various sizes 
around the origin, and so determines at least i — 1 layers of ordinary tiles around the origin. 
By choosing N large enough, we can get i to be arbitrarily large. Thus as N — > oo the 
points in the approximants determine tilings of larger and larger balls around the origin, so 
a point in the inverse limit gives a set of consistent directions for tiling increasing regions of 
M. d whose union is all of M. d . Such instructions are clearly in 1:1 correspondence with tilings 
of W*. Checking that the topology of X-r corresponds to the topology of the inverse limit 
(as a subset of the infinite product Y\ T n ) is a straightforward exercise that is left to the 
reader. □ 

Example 5.4. A short Cech cohomology computation. Consider a transition-regular fusion 
rule in one dimension, with two tile types a and b. Let P n (a) = P n _i(a)P n _i(fe)P n _i(6) 
and let P n (b) be a permutation of two P n _i(a)'s and three P n _!(6)'s, with the permutation 
depending on the level. As long as a permutation beginning in P n _i(a) occurs infinitely 
often and a permutation ending in P n _i(6) occurs infinitely often, this fusion rule forces 
the border. The approximant Y n consists of two intervals, one representing P n (a) and one 
representing P n (b), with all four endpoints identified to form a figure-8. The map f n wraps 
the P n+ i(a) circle around the P n (a) circle once and then around the P n (b) circle twice. It 
also wraps the P n+ i(b) circle around the P n (a) circle twice and around the P n (b) circle three 
times, in an order determined by the fusion rule at level n + 1. By Theorem 5.3, Xn is the 
inverse limit of these figure-8's under these maps. 

From a Cech cohomology standpoint we can see the figure-8 as the chain complex of 
each approximant, so that both Hi(T n ) and H l (T n ) are isomorphic to Z 2 . The first Cech 
cohomology of the inverse limit (and of Xn) is the direct limit of Z 2 under the pullback of 
the maps /„, which always come out to be (23) even though the order for the P n (6)'s varies. 
Since that matrix is invertible over Z, we see that H\X n ) = Z 2 . 

5.3. Anderson-Putnam collaring. Most fusion rules, like most substitution rules, do not 
force the border. However there is a simple trick, due to Anderson and Putnam in the setting 
of substitutions, for replacing an arbitrary fusion rule TZq with a hierarchical rule 1Z that 
forces the border, such that Xn is homeomorphic (and topologically conjugate) to Xn- We 
can then express our original tiling space Xn as the inverse limit of the Anderson-Putnam 
complexes of 71. 
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A collared tile to distance r, or an r-collared tile, is a tile together with a label that 
describes the types and relative positions of all of that tile's neighbors out to distance r. 
For instance, in a 1-dimensional tiling with patch abbaba, the three 6's are all different as 
collared tiles to distance 1, as one is preceded by an a and followed by a b, one is preceded by 
a b and followed by an a, and one is both preceded and followed by a's. Likewise, a collared 
n-supertile to distance r is an n-supertile, together with a label indicating the pattern of 
nearby n-supertiles out to distance r. 

Take an infinite increasing sequence of radii r Q < r± < ■ ■ ■ , tending to infinity. We take V n 
to be the set of collared n-supertiles to distance r n . Clearly, any complete tiling can be equally 
well-described in terms of (super)tiles or collared (super)tiles. However, by construction, 1Z 
forces the border, since if is greater than r n plus the diameter of the largest n-supertile, 
then a collared iV-supertile determines its surrounding r n -collared n-supertiles. 

Note that the label of a collared n-supertile contains information about all the neighboring 
n-supertiles out to distance r n , and in particular determines all of the neighboring (n — 
l)-supertiles out to distance r n _j. This means that each collared n-supertile is uniquely 
decomposed as a union of collared (n — l)-supertiles, and gives a well-defined map from 
to -X^T 1 . The hierarchical rule 1Z is a generalized fusion in the sense of Footnote 1, since the 
(n — l)-supertiles contained in an n-supertile do not determine the n-supertile. The collared 
n-supertiles have strictly more information than the collared (n — l)-supertiles, which is the 
whole point of collaring! 

If the shapes of the supertiles are not too wild, one can pick the r n 's to grow slowly 
compared to the size of the supertiles, so that collaring to distance r n just means specifying 
the nearest neighbors of the n-supertile, as is usually done for substitution tilings. However, 
for some fusion rules it is possible that knowing the n-supertiles containing the origin and 
the ones touching this supertile, for all n, will not determine the tiling of all of Mr. 

The process of collaring does have its drawbacks, as TZ may not have the same transition- 
regularity or even prototile-regularity properties as 71$. Collaring increases the number of 
tile types, and there is no reason to expect the increase to be the same at all levels. Indeed, 
even if the number of uncollared supertiles is uniformly bounded it is entirely possible that 
the number of collared n-supertiles will grow without bound as n — > oo. This happens, for 
instance, in Example 6.4. 

5.4. Barge-Diamond collaring. The idea behind Barge-Diamond collaring [6, 7] is to 
collar points rather than tiles. As before, pick an increasing sequence of radii r < r\ < ■ ■ ■ 
tending to infinity. Take a tiling T, and identify points x and y if [B ro ] T ~ x = [B ro ] y . 
That is, if the tiling looks the same around x and y to distance r (with x and y playing 
corresponding roles). Let r be the quotient space. To get Ti, identify points for which the 
corresponding tiling in X n agrees to distance r\. That is, points x and y for which all of 
the 1-supertiles that exist within a distance ri of x and y agree. Likewise, T n is M d modulo 
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identification of points x and y for which all of the n-supertiles within distance r n of x and 
y agree. 

As before, we have a map from X-ji to T n , taking a tiling to a description of how a ball 
of radius r n around the origin sits in one or more n-supertiles. Since r n — > oo as n — > oo, a 
point in the inverse limit is a consistent set of instructions for tiling all of M. , so the inverse 
limit is homeomorphic to Xn as long as the orbit closure of T is X-ji (which is always the 
case when X-ji is minimal). 

When the fusion is asymptotically self-similar or self-affine, one can take r n to be much 
smaller than the size of an n-supertile, but still to go to infinity. For 2-dimensional tilings, 
this means that there are three kinds of points. Most points are farther than r n from 
the nearest n-supertile boundary. These points are identified with corresponding points of 
other n-supertiles, without regard for the supertile's neighbors. Some points are within r n 
of one of the supertile's edges. These points are identified with corresponding points of 
other n-supertiles that have the same n-supertile neighbor across the specific edge. Finally, 
some points are within r n of two or more edges, and hence are close to a vertex. There 
is a stratification of T into points-near-vertices, points-near-edges, and interior points, and 
this stratification makes for much easier computations of tiling cohomology than Anderson- 
Putnam collaring. 



6. Direct product variations 

An easy way to make higher- dimensional substitution sequences is to take the direct 
product of two or more one-dimensional substitutions. To break the direct product structure, 
one can rearrange the substitution carefully so that at each stage the blocks still fit, creating 
what is called a direct product variation or DPV. Introduced as examples of combinatorial 
substitutions in [25], DPVs are quite flexible when viewed as examples of fusion rules. 



Example 6.1. The Fibonacci DPV. This simple example of a prototile- and transition- 
regular fusion rule in two dimensions is based on the Fibonacci substitution — > 01, 1 — > 0. 
We use it to illustrate almost all of the ideas and computations discussed for fusion rules. 
The prototile set consists of four unit-square tiles with label set {a, b, c, d} and so Vq = 



There is some flexibility with the geometry of the prototiles. They could be parallelograms or rectangles, 
and there are two vertical and two horizontal degrees of freedom for the lengths of the sides. 
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For the 1-supertiles we choose V\ 
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b 
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X 



where we list the super- 



tiles in the obvious order {Pi(a), Pi(b), -Pi(c), Pi(d)}. To make the 2-supertiles we concate- 
nate the 1-supertiles in combinatorially the same way: 
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In general we construct V n+ i from V n with exactly the same combinatorics as the rightmost 
version of the 2-supertiles shown above. It is not difficult to show that the sides of P n (a) 
and the long sides of P n {b) and P n (c) are the Fibonacci numbers f n +2, while the sides of the 
n-supertile of type d and the short sides of the h and c supertiles are the Fibonacci numbers 
f n+ i (using the convention that fo = and fi = 1). This means that at each stage, the 
supertiles fit together to form squares and rectangles with Fibonacci side lengths. 

Recognizability is straightforward and proceeds by induction. The P n+ i(a) supertiles are 
determined by the presence of a P n (d), each P n+1 (b) is determined by a P n (c) that is not in a 
P n +i( a ), each P n+1 (c) is determined by a P n {b) that is not in a P n+1 (a), and each remaining 
P n (a) is a P n+ i(d). 

(i i i i\ 

The Fibonacci DPV is transition-regular with M = M ra _i jn = } ? n o • This is a prim- 



itive matrix so M n ^ = M N ~ n is asymptotically rank 1; the dynamical system is uniquely 
ergodic. The Perron- Frobenius eigenvalue is 2 , where <fi is the golden mean; this number 
represents the asymptotic volume expansion of the supertiles from level to level. 

To compute the ergodic measure, it suffices to compute the frequencies of the n-supertiles 
and then use equation (3.4) of Theorem 3.4 to get the frequencies of arbitrary patches. The 
vectors p n are the volume-normalized directions of the asymptotic columns of M ny N = M N ~ n , 
and thus they are collinear with the right Perron-Frobenius eigenvector of M, which is 

. Since we have chosen unit square prototiles, the volumes of the n-supertiles are 

(<? 

/n+2> (/n+a/n+i)> (fn+2f n +i), and f% +1 respectively. We compute p n = (f) (2n+4) I J 

Next we turn to computing the topological and measure-theoretic spectrum. Technically 
we should take the induced fusion rule that composes two levels at once to get strong primi- 
tivity, then go up two more levels at a time to find all of the return vectors in V n . Fortunately, 
in this example there are always return vectors of the form (/„, 0) and (0, f n ) (see Figure 4). 
Any eigenvalue a = (04,02), topological or measure-theoretic, must have the property that 
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P n (a) 





Figure 4. The induced fusion for P n+2 (a). 

lim a-v n = (mod 1), and this means that lim a^fn — (mod 1) for k = 1, 2. Pisot's the- 
orem then implies that a G Z[0] x Z[0] 6 . For such a, the criterion for topological eigenvalues 
in Theorem 4.1 is satisfied because the convergence of rj n (a) to is exponential. In short, the 
computations for both topological and measure-theoretic eigenvalues are exactly the same 
as for the product of two Fibonacci tilings (built from unit length prototiles), and we have 
that both eigenvalue sets are Z[</>] x Z [</>]. In particular, any measurable eigenfunction can 
be chosen to be continuous. 

Next we compute the Cech cohomology of We use Barge-Diamond collaring [6, 7], 
picking the collaring radius r n to grow slowly with n. We stratify our tiling space into three 
pieces S C Si C S 2 = X*n- S 2 is the entire tiling space, Si is the set of tilings where 
the origin is within r n of the boundary of an n-supertile for every n, and So is the set of 
tilings where the origin is within a distance r n of two supertile edges, and hence is near a 
supertile corner. The cohomology of S is the cohomology of a CW complex with one cell 
for each possible pattern by which three or more high-order supertiles can meet at a vertex. 
There are 78 such patterns. The relative cohomology of the pair (Si, So) is computed from 
a CW complex containing 52 cells that describe the ways that two supertiles can meet along 
a common edge. The relative cohomology of the pair (S 2 , Si) is computed from the matrix 
M. The techniques for generating these cells and computing the cohomology are similar to 
those for substitution tilings, and yield 



^°(S ) 


= Z; 


^(So) = 


-0; 


^(So) = 


Z 42 


tf°(Si,S ) 


= 0; 


^(Si,So) 


= Z 4 ; 


# 2 (Si,S 
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if 2 (Si) = 


z 60 


^°(S 2 ,Si) 
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^(S 2 ,Si) 


= 0; 


i? 2 (S 2 ,Si 


) = z 4 


H°(X n ) 


= Z; 


H\X n ) = 


Z 4 ; 


H 2 (X n ) = 


= z 64 . 



The generators of H l (Xf>) = Z 4 are easily described. Pick a value of n > 4. Each edge 
of an n-supertile either has length f n+ \ or f n +2- The first generator counts the horizontal 
edges of the first type, the second generator counts the horizontal edges of the second type, 

6 The absence of the V5 that is present in Example 4.4 is due to the integer size of the prototiles. 



42 



NATALIE PRIEBE FRANK AND LORENZO SADUN 



and the third and fourth generators similarly count vertical edges. The boundaries of two 
supertiles may overlap on intervals of size /„ = f n+2 — f n +i or / n _j = 2f n+1 — f n+ 2- The first 
(or third) generator assigns the numbers —1 and 2 to these partial edges, while the second 
(or fourth) assigns the numbers 1 and —1, as these are the coefficients of f n+ i and /„+2- 
Picking different values of n gives different generators, but the group they generate is the 
same. 

Deformations of a tiling, by changing the shape and size of (possibly collared) tiles, are 
parametrized up to mutual local derivability by iy 1 (X^, W 1 ) [16]. For the Fibonacci DPV, 
H 1 is the same as for the product of two 1-dimensional Fibonacci tiling spaces, and the 
deformations are the same. Thus, any deformation of the sizes and shapes yields a tiling 
space that is topologically conjugate to a linear transformation of 1R 2 applied to the original 
tiling space. In particular, a self-similar version of the DPV, in which the a, b, c, and d tiles 
have dimensions 0x0, 0x1, 1x0 and 1 x 1, is topologically conjugate to a DPV where 
all tiles are congruent squares (of side \/5/0). 

The difference between the Fibonacci DPV and the product of two 1-dimensional Fibonacci 
tilings is seen in the second Cech cohomology, where that of the DPV has rank 64 and that of 
the product has rank 4. The rank of the top cohomology is closely related to the independent 
appearance of patterns in the tiling, via the following theorem: 

Theorem 6.2 ([ ]). If the rank of H d of a d- dimensional tiling space is k, then there exist 
k patterns Pi, . . . , Pk, such that for any patch P there exist rational numbers Cj, . . . , c& and 
Cp such that, for any region R in any tiling T , 

k 

#(P %n R) = 5> t #(P, in R) + e(P,R), 

i=l 

where the error term e(P, R) is computable from the patterns on the boundary of R, and is 
bounded by cp times the (d — 1) -volume of the boundary of R. 

We call P\, . . . , Pk control patches. For the product of two 1-dimensional Fibonacci tilings, 
we can take our control patches to be the four basic tiles. For the Fibonacci DPV, however, 
there are 60 additional control patches. They can be chosen from the generators of H 2 (E,q) 
and H 2 (Ei,E ). That is, we have 9 control patches that involve supertiles meeting along 
horizontal edges, 9 that involve supertiles meeting along vertical edges, and 42 that involve 
three or four supertiles meeting at a vertex. 

Example 6.3. A scrambled Fibonacci DPV. We can construct a scrambled version of the 
Fibonacci DPV in much the same way as the 1-dimensional scrambled Fibonacci tiling of 
Example 4.4. We pick an increasing sequence N(n) and induce on this sequence to get an 
accelerated scrambled Fibonacci rule A. We then introduce an exceptional supertile £ n (e) 
at each odd level, whose population in terms of (n — l)-supertiles is the same as A n (d), but 
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rearranged so that all of the S n ^i(a) tiles appear in the lower left corner, all the £ n _x(6) 
appear in the lower right, all the £ n _i(c) appear in the upper left, and all the £ n _x(d) appear 
in the upper right. On even levels, the n-supertiles are built from the (n — l)-supertiles 
exactly as for the accelerated DPV, only with one S n -i(d) in each ra-supertile replaced by 
an £ n _i(e). Finally, we induce on even levels to obtain a prototile-regular fusion S. 

As before, if we choose the sequence N(n) to grow sufficiently fast, and if we give the 
prototiles the same shape as the asymptotic supertiles, with the a, b, c, d prototiles having 
dimensions 0x0, 0x1, 1x0 and lxl, then the scrambled Fibonacci DPV space is 
topologically weakly mixing but has pure point measurable spectrum, being measurably 
conjugate to the unscrambled Fibonacci DPV. However, if we choose the prototiles to be 
unit squares, then every a G Z x Z is manifestly a topological eigenvalue. 

This discrepancy means that the deformation theory for the scrambled DPV is not the 
same as for the unscrambled DPV. Either the first cohomologies are different, or, more likely, 
the cohomologies are isomorphic but the two tilings have different "asymptotically negligible" 
[16] subspaces of H 1 (X, M. 2 ) that describe deformations that are topological conjugacies but 
that are not mutually locally derivable from the original. As for H 2 , the rank must be at 
least 64, since the control patches for the Fibonacci DPV are still present in the scrambled 
DPV. 

This example suggests two directions for future work. One is to understand deformation 
theory better, and in particular the role of the asymptotically negligible classes. These are 
well-understood for substitution tilings, but not for fusions. Another is to develop new tech- 
niques for computing tiling cohomology for spaces that do not come from substitutions. The 
Anderson-Putnam and Barge-Diamond complexes were defined in Section 5 for all tilings, 
but almost every existing method for studying these complexes relies on an underlying sub- 
stitution. 

Example 6.4. A non-Pisot DPV. We base this DPV on the one-dimensional substitution 
a — > abbb, b — > a, which despite its apparent similarity to the Fibonacci DPV exhibits 
significantly different dynamical behavior. The prototile set is the same as for the Fibonacci 
DPV, and again we choose the fusion to be both prototile- and transition-regular. This time 
we choose our fusion rule at each stage to be given by 
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Recognizability is easily established, almost exactly as with the Fibonacci DPV. 

The transition matrix has ((1 + V4~3)/2) 2 as its largest eigenvalue, which is not a Pisot 
number. The side lengths of the supertiles grow as solutions to the recursion Z n+1 = l n + 
3s n ,s n+ i = l n , which are nontrivial linear combinations of ((1 ± vl3)/2) r \ Since both of 
those numbers are greater than one in modulus, the side lengths are not well-approximated 
by powers of the positive eigenvalue. The effect this has on the system is profound. It means 
that the combinatorics of the fusion tilings are exceptionally complicated, in that the number 
of ways that n-supertiles can be adjacent to one another grows without bound as n — > oo. 

This increasing complexity with scale shows up in the topology of Both H 1 and 

H 2 are infinitely generated, the first indicating that there are infinitely many "interesting" 
deformations of size and shape, and the second indicating that there are infinitely many 
control patterns. Unlike the Fibonacci DPV, changes in tile sizes, while preserving the 
fusion rule, can change the dynamics and in fact the topology of the tiling space. If we were 
to choose irrationally related side lengths for our prototiles, then the resulting tiling would 
not have finite local complexity [26]. 

Examples 6.1 and 6.4 lead us to a discussion of the combinatorial and geometric behavior 
of supertiles as n — > oo. In some cases one or the other will approach a limit as n — > oo. 
Consider a prototile-regular fusion rule, and suppose that there is some invertible linear 
map L : IR d — > IR d such that lim L~ n P n (p) exists for each prototile type. If in addition 
the combinatorics of how the (n — l)-supertiles lie inside their n-supertiles stabilizes for 
large values of n, then we call the fusion rule asymptotically self-affine (or -similar if L is a 
similarity). This means that there is a self-affine tiling that is related to the fusion tiling. 
The precise nature of the relationship varies, and no general theorems about it are known 
to the authors at this time. Both of the previous examples are asymptotically self-similar, 
with the limiting prototile sets having edge lengths in the ratios <fi : 1 in the Fibonacci case 
and (1 + V4~3) : 2 in the non-Pisot case. 

A fusion tiling may have finite local complexity in the usual sense while failing to be locally 
finite in an asymptotic sense. We call a fusion rule asymptotically FLC if there is a constant 
B such that each pair of n-supertiles can form at most B connected two-supertile patches. 
Example 6.1 is asymptotically FLC, but Example 6.4 is not. If an asymptotically self-affine 
tiling is not asymptotically FLC, then the self-affine tiling obtained from the limiting shapes 
will have infinite local complexity. 
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